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Several  Improvements  In  previous  treatments  of  the 

relaxation  theory  of  spectral  line  broadening  In  a  plasma 

are  presented  in  this  dissertation i 

1.  To  second  order  in  the  coupling  constant «  the  theory 
is  extended  to  include  ion- field  effects  in  the  atomic 
Boltzmann  factors  and  the  electron  collision  operators. 

2.  The  theory  is  extended  to  include  full  Coialomb 
interaction  betveen  the  radiating  atoms  and  the  perturbing 
electrons. 

3.  The  effect  of  using  a  strong  collision  cut-off  on 
the  theoretical  computations  is  examined. 

4.  The  effects  of  electron  symmetry  are  Introduced 
through  an  effective  Boltzmann  distribution. 

5.  A  static  wing  formula  is  developed. 

Computations  are  presented  for  the  Lyman-alpha  transition 

vll 


in  hydrogen  and  ionized  helium*  For  frequency  separations 
from  the  line  center  greater  than  10  times  the  ion  plasma 
frequency,  the  calculated  profiles  for  the  Lyman-alpha  line 
agree  quite  veil  with  the  experimental  data*  Hovever,  there 
are  noticeable  differences  in  the  theoretical  results  in  the 
line  center  I  these  discrepancies  are  examined* 

From  the  behavior  of  the  Lyman-alpha  profile,  when  the 
full  Coulooib  interaction  is  used,  it  is  concluded  that 
higher  order  terms  in  the  coupling  constant  are  required  to 
properly  describe  the  Intensity  in  the  line  center*  An 
explicit  expression  is  derived  for  a  typical  higher  order 
term  in  the  perturbation  expansion  and  an  upper  bound  is 
obtained  for  this  term.  It  is  shown  that  the  perturbation 
expansion  does  indeed  break  down  for  small  frequency 
separations  corresponding  to  the  line  center* 

The  calculated  results  using  the  static  ving  formula 
are  compared  vith  the  experimental  data  for  frequency 
separations  greater  than  the  electron  plasma  frequency. 
This  theory  provides  an  Inexpensive  means  for  computing  wing 
profiles  which  may  be  useful  to  experimentalists  and  astro- 
physicists* 
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CHAPTER  I 
CAUSES  OF  SPECTRAL  LINE  BROADENING  IN  PLASMAS 

(I.l)  INTRODUCTION 

The  discrete  spectrum  of  radiation  emitted  from 
isolated  atoms  and  molecules  is  broadened  considerably  vhen 
the  radiators  are  placed  in  a  plasma*  The  resulting 
intensity  distribution  as  a  function  of  frequency  is  called 
a  line  shape*   The  major  caiuses  of  line  broadening  are 
discussed  in  Chapter  I*  The  detailed  behavior  of  the  line 
shape  under  different  approximations  is  investigated  in 
sxibsequent  chapters. 

The  broadening  mechanisms  can  be  grouped  into  three 
major  categories i  the  effects  of  charged-particle  fields  on 
the  radiation*  collisions  between  the  radiator  and  various 
neutral  constituents  in  the  plasma*  and  Doppler  effects* 
In  this  chapter,  the  Lyman-alpha  line  of  hydrogen  vill  be 
considered  in  order  to  illustrate  the  region  of  the  line 
shape  vhere  the  different  broadening  effects  are  important} 
unless  othervlse  stated,  all  numerical  results  pertain  to 
a  hydrogen  plasma  at  a  temperature  of  20,400  degrees  Kelvin 
and  an  electron  density  of  3.6  x  lO^^/cm^.  However,  most  of 
the  analysis  will  hold  for  temperatures  and  densities  which 
may  vary  by  factors  as  large  as  4  for  the  temperatures  and  10 
for  the  density!  exceptions  will  be  noted* 
.  „        1 


(1.2)   TIME  OF  INTEREST 
To  examine  the  effects  of  the  various  broadening  mech- 
anisms on  the  line  shape.  It  Is  useful  to  define  a  time  of 

2 
interest.   The  line  shape  can  be  expressed  as  a  Fourier 

transform  of  an  auto-correlation  function*  ^(X.)»  '  ' 

vhere  $(je)  «  "1^^ JL»  3f*>^  andAO;  is  the  frequency 
separation  from  the  line  center*  For  times  such  that^^Z^l, 
contributions  to  the  Integral  are  small  due  to  oscillations 
of  the  exponential  function.  Also,  for  cases  investigated^ 

in  this  dissertation,  ^(Z)   can  be  shown  to  be  small  for 

5 
large  times.   Thus,  contributions  to  the  integral  are 

significant  only  for  times  less  than  1/^kW  .  Thus,  t.  "  1/Ck.^ 

defines  the  time  of  interest. 

Experimental  profiles  for  the  Lyman-alpha  line  for 

wavelengths  less  than  the  half  width,  /iX^  ■»  'IS  £,  are 

6  7 

not  well  defined.  '   Using  this  value  of  ^X^^  the 
corresponding  time  of  interest  can  be  calculated. 

Changes  in  physical  quantities  such  as  the  ion-field,  which 
are  appreciable  only  for  times  greater  than  '^^u/^  have  small 
influence  on  the  profile  for  ^>><^>^j.^.  Hence,  they  will 


3 
be  neglected  In  this  dissertation i  for  example,  the  static 
ion  approximation  vill  be  employed. 


(1.3)  DOPPLER  BROADENING 

Doppler  broadening  Is  caused  by  motion  of  the  radiator 
during  the  time  of  Interest,^  Neglecting  the  effects  of 
collisions  and  the  ion-field  on  the  motion  of  the  radiator, 
the  Doppler  profile  can  be  written  as. 


.'«)44?)^-?^)' 


(1.1) 


vhere  M  is  the  mass  of  the  radiator,  c  the  velocity  of  light, 
J^  the  Bolt zmann  constant,  and  Ap  the  Debye  length.  For 
cases  considered  in  this  dissertation,  the  Doppler  broaden- 
ing vlll  be  assumed  statistically  Independent  from  other 
broadening  effects.  Under  this  assumption  the  line  profile 

o 

can  be  represented  as  a  convolution  of  the  Doppler  profile 
and  the  profile  resulting  from  all  other  broadening  mechanisms, 

I^^,(^.^)=  5*I/AA')I^{-iA-AA')d(AX)»         a.2) 

tjhere  lOiX  )  is  the  profile  arising  from  all  other  effects. 

For  long  times  of  Interest,  tdiich  correspond  to  the 
line  center,  the  motional  effects  of  the  radiator  will 
become  more  pronounced.  Thus  In  most  cases,  Doppler  effects 
are  significant  only  in  the  line  center.  This  fact  is 
Illustrated  in  reference  9. 


5 

The  Doppler  profile  given  In  equation  (l.l)  also 
neglects  the  effects  of  electron  collisions  on  the  motion 
of  the  radiator. ^°'^^  This  is  a  valid  approximation  since 
the  momentum  transferred  to  the  radiator  during  an  electron 
collision  is  small  due  to  the  large  difference  in  the 
electron  and  atom  masses* 

In  this  dissertation,  attention  vill  be  focused  on 
theoretical  treatments  of  K^X  )• 


(1.4)  PERTURBING  IONS  AND  ELECTRONS 

Almost  all  treatments  of  line  broadening  to  date  use 

1  12  13 
the  static  Ion  approximation*  '   '    In  this  approximation 

the  broadening  effects  of  ion  motion  or  ion  dyneunics  are 

neglected*  In  reference  9  it  is  shown  that  the  relative 

motion  of  the  ions  is  small  for  times  corresponding  to 

wavelengths  greater  than  the  plasma  frequency  for  the  ions* 

For  the  cases  investigated  in  this  dissertation*  the  plasma 

frequency  for  the  ions  is  smaller  than  the  half  width* 

Hence,  the  effects  of  ion  dynamics  can  be  neglected  for 

wavelengths  greater  than  ^X^^jj   • 

The  criterion  used  in  determining  the  region  of  the  line 
shape  in  which  electron  motion  is  Important  is  the  same  as 
that  used  >^en  ion  dynzunlcs  are  considered*  The  wavelength 
corresponding  to  the  electron  plasma  frequency  is  2  &*  Thus 
electron  dynamics  will  have  relatively  less  importance  for 
wavelengths  greater  than  2  X*  As  will  be  shown  in  Chapters 
IV  and  V»  a  complete  static  treatment  for  both  the  electrons 
and  ions  yields  a  theoretical  profile  that  is  in  substantial 
agreement  with  the  experimental  data  in  the  line  wings* 

The  effects  of  electron  shielding  on  the  ion  mlcrofleld 
distribution  function  must  also  be  considered*  A  comparison 
of  the  relative  motion  of  the  electrons  to  that  of  the  ions 
during  the  time  of  interest  yields 


(^•^A*^Xi.»cr«i>»/ 


[hhj^lJ^Js^ 


4^^ 


43. 


•LO^ 


Since  electron  motion  during  the  time  of  interest  is  large, 
the  electrons  may  be  considered  to  effectively  shield  the 
ion-ion  interaction.  Thus  the  appropriate  ion  microfield  is 
the  lov^frequency  microfield.^ 

The  same  argument  can  not  be  applied  to  a  treatment  of 
the  electron  microfield.  Because  of  the  lov  mobility  of 
the  ionSf  the  electron-electron  interaction  is  not  fully 
shielded.  To  explore  the  effects  of  shielding  in  the  elec- 
tron microfield,  the  static  wing  profile  calculated  using  a 
low-frequency  electron  microfield  is  compared  to  that  obtained 
using  the  high-frequency  electron  microfield. ^^  These  tvo 
cases  should  provide  bounds  for  the  eflbct  of  ion  shielding 
on  the  electron  microfield. 

In  the  line  center,  0.0  -  5.0  X,  for  the  temperature 
and  density  under  investigation,  the  interaction  between  the 
radiating  atom  and  the  perturbing  electrons  is  treated  as  a 
collision  process  (i.e.,  no  static  electron  field  broadening 
is  Included) .  One  approximation  often  used  in  describing 
these  collisions  is  to  assume  that  the  electrons  travel  in  a 

13   16 

classical  path.      '  in  many  of  these  treatments  logarith- 

mic divergences  are  encountered  when  the  dipole  approximation 
is  used,   and  only  terms  to  second  order  in  the  coupling 
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17  18 
constant  are  retained.   '    These  divergences  are  removed 

by  using  minimum  and  maximum  impact  parameters. ^^'^°'^^'^^ 

The  maximum  impact  parameter  is  assumed  to  be  the 

Debye  length  or  a  constant  close  to  unity  times  the  Debye 

length.  Several  different  criteria  are  used  to  determine 

the  minimum  impact  parameter  cut-off  or  strong  collision 

cut-off.  For  the  Lyraan-alpha  transition,  the  minimum 

impact  parameter  is  often  in  the  region  in  which  the  dipole 

approximation  for  the  Interaction  between  the  radiator 

and  perturbing  electrons  is  inadequate.  By  using  the  quantum 

mechanical  treatment  developed  in  the  relaxation  theory, 

employing  either  the  dipole  approximation  or  the  full  Coulomb 

interaction,  there  is  no  mathematical  need  for  a  strong 

collision  cut-off  in  the  second  order  term  for  frequencies 

greater  than  the  half  width.  However,  in  order  to  obtain 

realistic  results  a  strong  collision  cut-off  is  required, ^^ 


(I ,5)  INTERACTIONS  BETWEEN  THE  RADIATORS 

Collisions  betveen  the  neutral  radiators  are  neglected 
in  roost  plasma  line  broadening  treatments.  To  examine  the 
validity  of  this  approximation*  the  nature  of  the  radiator- 
radiator  interaction  must  be  examined.  One  can  compute  the 
average  density  of  neutral  atoms  from  a  knowledge  of  the 
percent  ionization.  The  percent  ionization  can  be  determined 
from  recent  computations  by  Barker.    For  an  electron 
density  of  1.0  x  10  /cm  and  a  temperature  of  20,000  degrees 
Kelvin,  the  percent  ionization  is  approximately  75  percent 
corresponding  to  an  average  distance  between  radiators  of 
approximately  1.0  x  10   cm.  During  the  time  of  interest 
the  ratio  of  the  average  motion  of  the  radiator  to  the 
average  separation  is  less  than  1.  Also,  since  the  radiator- 
radiator  interaction  potential  is  short  range,  typically  a 

6 
1/r  dependency,  the  effective  mean  free  path  will  be  much 

25 
larger  than  the  distance  between  radiators.    Thus,  during 

the  time  of  interest  there  will  be  few  radiator-radiator 

collisions.  The  broadening  due  to  these  collisions  will  be 

neglected  compared  to  the  effect  of  the  electron-radiator 

collisions  which  interact  through  long  range  potentials. 


10 


(1.6)  THE  NO-QUENCHING  APPROXIMATION 

In  the  context  of  line  broadening  theories  the  no- 
quenching  approximation  neglects  radiationless  transitions 
betveen  the  initial  atomic  states  and  all  other  atomic 
states  having  a  different  principal  quantum  number.  This 
is  a  good  approximation  betveen  states  tihich  are  well 
separated  in  energy  such  as  the  lower  transitions  in  the 
first  few  series  in  hydrogen  and  helium.  However,  quenching 
can  become  appreciable  for  higher  members  in  a  series  for 
which  the  transition  probability  to  states  of  different 
principal  quantum  numbers  is  large.  Including  quenching 
effects  in  the  lover  member  transitions  is  computationally 
feasible  I  unfortunately,  this  is  not  the  case  for  higher 
transitions  because  integrations  over  the  inverse  of  large 
matrices  are  required.  Quenching  has  been  examined  for  the 
Lyman-alpha  transition  by  using  a  matrix  perturbation 
expansion  approach*   It  Is  shown  that  the  effect  is  small 
(see  Chapter  V).  The  approach  used  could  be  extended  to 
the  higher  series  members.  However,  at  this  time  it  is 
felt  that  such  an  extension  would  be  premature  since  other 
approximations  have  a  more  significant  effect  on  the  line 
profile  and  should  be  removed  first. 


CHAPTER  II 
THE  RELAXATION  THEORY 

(II.l)  INTRODUCTION 

9  12 
The  relaxation  theory  '   of  spectral  line  broadening 

in  plaamas  la  reformulated  In  this  chapter  to  include  the 

following  I   (1)  the  full  Coulomb  Interaction  betveen 

radiating  atoms  and  perturbing  electronst  (2)  lon-fleld 

splitting  of  Initial  atomic  levels*  both  In  the  Boltzmann 

factors  and  In  the  electron  collision  operators. 

In  the  relaxation  theory  the  profile,  l(^)»  of  a 

broadened  spectral  line  Is  given  by  the  Fourier  transform 
^^9,12,26 


Hence » 


(2.1a) 


(2.1b) 


vhere  Tr  represents  a  trace  over  all  states  of  the  plasma 
system  vhlch  consists  of  radiators*  electrons*  and  Ions. 
T(t)  Is  the  time  development  operator  and  ^  Is  the  density 
matrix  of  the  system.  This  dissertation  Is  concerned  with 
atomic  dlpole  radiation  from  partially  Ionized  plasmasi 
hence*  Ji   Is  restricted  to  be  the  atomic  dlpole  operator  of 
the  radiator. 

11 


12 
Uhder  the  conditions  outlined  in  Chapter  I,  the  radia- 
tors can  be  asstoned  isolated  from  one  another.  Also  all  the 
radiators  vill  have  the  same  spectrum  of  states  and  are 
independent  of  one  another  as  shovn  in  Section  (I* 5).  Hence, 
to  obtain  a  useful  line  profile  expression,  it  is  sufficient 
to  consider  the  ensemble  average  over  all  states  accessible 
to  a  single  radiator*  The  corresponding  Hamiltonian  for  a 
radiator  immersed  in  a  gas  of  charged  particles  is 

H=  H.  +  H*  +  K,.  +  V.v  +\4..  t-  M^^  ,2,2, 

t.There  subscripts  a,  e,  i  represent  atoms,  electrons,  and  ions, 
respectively.  V  represents  the  total  potential  energy  opera- 
tor, and  K  the  total  kinetic  energy  operator.  H^  is  the 
Hamiltonian  for  the  electrons  defined  by 

9  12 
Previous  treatments  '   of  the  relaxation  theory  used  an 

unperturbed  atomic  Hamiltonian  of  the  form  Kn  •«■  V^^.   In  this 

dissertation,  the  effect  of  the  ion-field  on  the  atomic 

states  is  included  and,  hence,  the  unperturbed  atomic 

Hamiltonian  in  the  following  derivation  is  taken  to  be 

vhere  e  is  the  electron  charge,  R  the  position  operator  for 
the  atomic  electron,  and  £  .  is  the  static  ion-field. 

To  at  least  qualitatively  account  for  electron  shielding 
(see  Chapter  I,  Section  1.4)  the  operator,  V. .  +   V  4>  is 


13 
replaced  by  an  effective  potential »  V^ff*  vhlch  depends  only 

on  Ion  coordinates.  Debye  HUckel  theory   for  Ions  In  a 

neutralizing  electron  background  Is  used  to  obtain  V  ^ .. 
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(II. 2)  THE  DENSITY  MATRIX 


The  canonical  density  matrix  for  a  ■ystem  vhose 

Hamlltonlan  Is  given  by  equation  (2.2),  vlth  V^^  +  V^^ 

28 
replaced  by  V  -^f  Is 

vhere  -^le  Boltzmann*s  constant,  T  Is  the  temperature  and 
c  Is  a  normalization  constant  obtained  by  requiring  that 

Making  the  static  Ion  approximation  Is  equivalent  to 
assuming  that  the  following  commutator  Is  negligibly  small 

f;+V>  «]   =  "•  (2.4) 

since  K.  +  ^^ff  commutes  vlth  Itself  and  commutator  rela- 
tions are  distributive*  equation  (2.4)  Implies  that 

[H.+  /J,t\'«,  1^+- Vl   =  °-  «.5) 

Using  equation  (2.5),  ve  may  write  the  density  matrix 
In  the  following  convenient  form 


13 
The  operators  V_^  and  H  will  not  In  general  commute  because 
V   contains  the  atomic  electron  position  operator  and  H 
contains  the  atomic  electron  momentum  operator*  Ho^weverf  the 
exponential  factor  containing  V^^  may  be  expressed  as 


vhere  V   is  an  effective  potential  t^hich  is  defined  to 

depend  only  on  the  electron  coordinates i  an  average  over 

atomic  electron  coordinates  is  performed  to  obtain  vl^.  it 

is  assumed  that  the  term  (V_.  >  V*  _)  is  small  compared  to 

«"    ae 

both  V   and  ^,     Hence »  we  neglect  it  in  a  first  approxi- 
mation allowing  equation  (2.7)  to  be  factored* 

^^        *  ^  6  -"♦    (2.8) 


The  entire  density  matrix  can  nov  be  written  as  a  product 
of  the  density  matrices  for  the  three  subsystems  •••  atoms* 
electrons*  and  ions 


vhere 


•«^    ^^        *  (2.9a) 

fij  ^  ^  *  (2.9b) 

C        =^   normalization  constant. 
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We  next  consider  further  approximations  to  f^     .  It 
Is  often  stated  or  asstuaed  that  veak  electron  collisions 
account  for  the  major  contribution  to  the  electron  broaden- 
ing. However,  it  is  not  usually  pointed  out  that  a  primary 
reason  why  a  perturbing  electron  is  relatively  unlikely  to 
get  close  enough  to  a  radiator  to  cause  a  strong  collision 
is  due  to  electron  synunetry  effects*  For  the  purpose  of 
estimating  symmetry  effects  between  the  atomic  electron 
and  perturbing  electrons,  an  atomic  model  is  adopted  in 
which  the  atomic  electron  is  superimposed  on  the  nucleus, 
with  this  rough  approximation,  V*  ^  can  be  written  as  a  sum 
of  the  classical  Coulomb  interaction,  V  , ,  plus  a  term  V 

ex  g 

arising  from  symmetry  effects.  Since  the  classical  poten- 
tial produced  by  the  above  model  is  zero,  v"  reduces  to  V  . 

ae  s 

The  expression  for  the  density  matrix  f>    is  factored 
yielding 


p^  =  jt^ir  n^/Jkr)j>4f^(''  %/hT ) 


(2.10) 


where  the  commutators  between  H^  and  V^  are  neglected.  Again, 
in  the  spirit  of  the  previous  approximations,  it  is  assumed 
that  this  commutator  will  be  small  compared  to  ^T.  For 
temperatures  under  consideration  in  this  dissertation,  the 
form  for  ^  given  in  equation  (2.10)  should  be  adequate 
when  estimating  the  qualitative  effects  of  electron  symmetry. 

The  potential,  v^,  is  determined  from  the  radial  dis- 
tribution functions  for  electron-proton  and  an  electron- 

24  29 

electron  system,  as  computed  by  Barker.   '    The  details  for 
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deriving  the  radial  distribution  functions  and  the  corre- 
sponding effective  potentials  are  discussed  In  Appendix  A. 

A  plot  of  exp(-Vg/^T)  versus  the  distance  betveen  the 
perturbing  electron  and  radiator  Is  given  In  Figure  1.  The 
rapid  decrease  In  magnitude  for  small  distances  Implies  a 
decreased  probability  for  strong  collisions.  The  resulting 
effect  on  the  line  profile  can  be  represented  by  an  effective 
strong  collision  cut-off.  The  value  of  the  cut-off  obtained 
in  this  manner  is  indicated  in  Figure  1.   It  is  larger  than 
cut-offs  used  in  other  treatments  of  line  broadening.  Thus, 
for  this  model  electron  symmetry  effects  tend  to  decrease 
the  broadening  caused  by  the  electron  collisions. 
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(II. 3)  THE  TIME  DEVELOPMENT  OPERATOR,  T(t) 

The  time  development  operator  appearing  in  ^(T>)   ia 
defined  in  terms  of  the  Hamiltonlan  given  in  equation  (2*2) 
with  Vg£  +  ^ii  J^eplaced  by  V   . 

(2.11) 

vhere  the  commutator  relation  given  in  equation  (2.5)  has 
been  used.  Using  the  commutator  relations  given  in  equa- 
tions (2.4)  and  (2.5)  one  can  show  that 

Tif^ftT(f>  =/?-^(-'  i*WJi)fe  3  ^(/wV/4)       ^^  .^^j 


where 


(2.13) 


T    ( ^)  f^  f^  ^  T  i  ^)  is  evaluated  in  Section   (II. 5)   using 
the  Zwanzig  projection  operator  technique. 
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(11,4)  THE  AVERAGE  OVER  THE  ION  MICROFIELD 

Using  the  properties  of  the  Dirac  delta  fimction*  ve 
can  vrite  $(;(*)  as 


vhere  QUE.)  "=■    7^-  (f/  ^^^'^/'))   is  the  ion  microfield 
distribution  function.    Using  the  above  form  for  3  (f) 

in  the  expression  for  the  line  profile  given  in  equation 

9  12 
(2.1b),  one  can  show  that 

Ko))  =  if 'l?uuue  5"<V(^)  f^  ^^L^(J<o:r)Jju^  {t^-iW:f/j,) 

o 

It  should  be  emphasized  that  here,  unlike  previous  treatments 
of  the  relaxation  theory,  the  ion-field  is  retained  in  H*. 
Including  the  ion-field  introduces  field  shifted  frequencies 
in  the  effective  Liouville  operator  occurring  in  iT'i  <w ,  £  ) . 
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(II. 5)  APPLICATION  f)T   THE  ZWANZIG  PROJECTION  OPERATOR  TECHNIQUE 

We  must  nov  reduce  the  line  profile  expression  to  a 
tractable  form*  A  useful  technique  to  effect  this  reduction 

employs  the  projection  operator  procedure  developed  by 

4 
Zwan2sigi  this  procedure  involves  defining  a  projection 

operator «  P»  that  projects  out  the  relevant  part  of  an 

operator »  D(X),  that  satisfies  an  equation  of  the  form 

tihere  j_^  is  the  Llouville  operator.  The  projection  operator 
allows  one  to  remove  the  irrelevant  part  of  Ddf")  at  an  early 
stage  in  the  derivation*  thus  often  simplifying  the  analysis. 
For  example,  if  D.  is  the  relevant  part  of  D  and  D2  is  the 
irrelevant  part,  one  has 

■^(t^  =    P  D(^),  (2.17) 

D^U)=:  (l-P^b'<^).  (2.18) 

Applying  the  projection  operator  procedure  results  in  an 
integro-dif f erential  equation  for  D  (  f)  • 

It  is  often  unclear  exactly  vhat  part  of  an  operator  is 
relevamt  particularly  during  the  earlier  steps  of  a  deriva- 
tion. A  simple  example  of  a  relevant  part  of  an  operator 
might  be  the  diagonal  elements  of  a  matrix  operator.  A 
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projection  operator  vhlch  selects  only  the  diagonal  elements 

of  a  matrix  operator  vould  be  appropriate  -when  the  non- 
diagonal  matrix  elements  are  negligibly  small  or  vhen  the 
non-diagonal  elements  in  no  way  affect  the  final  result.  An 
example  of  the  latter  situation  is  vhen  an  average  of  an 
operator.  A,  is  taken  in  the  following  manners 

The  non-diagonal  elements  of  (^A)  in  no  way  effect  the 
average  value  of  the  operator  A. 

In  the  case  of  line  broadening,  one  is  concerned  with 
the  behavior  of  the  radiator,  either  an  ion  or  an  atom  in 
the  cases  considered  in  this  dissertation.  Hence,  the 
atomic  subsystem  can  be  considered  the  relevant  subsystem. 
One  performs  an  average  over  the  broadening  effects  of  the 
ions  and  electrons.  The  ion  broadening  effects  are  taken 
into  account  via  the  microfield  average.  A  suitable  proce- 
dure for  simplifying  the  electron  average  involves  a  projec- 

I       9 
tion  operator  of  the  following  form 

P-  Ci'^\  (2.19) 

where   ^  is  defined  in  equation  (2.10).   It  should  be  noted 

2 

that  P  is  indempotent,  i.e.,?^?. 

The  projection  operator  technique  Is  now  applied  to  the 
expression  appearing  in  equation  (2.12).  we  define  as  D('£  ) 

"D^:<^^"  Tjt)  fiT(^)  (2.20) 
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which  is  the  Helsenberg  representation  for^Jl^^ 

The  procedure  for  finding  an  expression  for  D^(  t")   is 
to  apply  the  operator  P  and  (1-P)  to  equation  (2.16)  and 
then  eliminate  the  operator  D-(  ^  from  the  resulting 
equations.  In  this  manner  one  obtains  an  integro-differential 
equation 

€>6 


-'i^  Wk-^<'-p>lM)0-P)1.^,<;^-^;)  '^^. 


For  most  problems  of  interest  the  operator  02(0) 


(2.21) 


occurring  in  equation  (2.21)  is  zero.  This  is  shoim  for 

9 
the  projection  operator  P  defined  in  equation  (2.19) 


^ 


(2.22) 


Thus  equation  (2.21)  becomes 

c-Jr:^  =  pl[d.(;^V  c'/^  ^^^(-/<uO-P>M) 
,^,-P)L'S,(^-^>A^]^ 

-*  (2.23) 

The  following  integral  must  be  determined  to  evaluate 
the  line  shape  expression 1 
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Ob  -A 


(2.24) 


■where 


An  equation  for  ^^(<y)  can  be  obtained  by  taking  the 
Fourier  transform  of  equation  (2.23).  Ualng  the  procedure 
outlined  In  reference  9  one  can  show  that 


P 

(2.25) 


<     >    Is  defined  by   ^A^    ■      T;»^  (A  f^   )  vhere  A  Is  an 

arbitrary  operator.     T (45 )   satisfies  the  Llppmann-Schwlnger 
equatlon^^  _^i 

(2.26) 

Substituting  the  above  result  Into  the  expression  for 
the  line  profile  given  In  equation  (2.15)  yields 
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j<^^,  Tf'x«,  Sw) \  0  •  h-«-i,^«)  ^  )  «•'£ 

*  (2.27) 

Inhere  the  resolvent  operator  K(u))  has  been  introduced  i 

|t'(^>  -   (<a>-  <T<**)>).  (2.28) 

The  Llouville  operator*  La  encountered  in  the  above 
derivation  is  defined  in  terms  of  the  Hamiltonian*  H', 
specified  in  equation  (2.13)  and  the  identity  operator  Ii 

L-   fl»3:  -X<«H^  (2.29) 

^     stands  for  a  direct  product.  L  can  be  rewritten  as  a 
sum  of  tvo  operators  Lq  and  L^^   defined  by, 

1^     **  '    '-^^  (2.30) 


L=(^^^^*^tMe)»>^-  x»?^/^--^^*.'^  %y 


L^     -        ""ajf  ^  *<  .  (2.31) 

Unlike  the  derivation  in  reference  9,  the  ion-field  depen- 
dence is  retained  in  L^.   If  equation  (2.26)  for  T(a>)  and 
the  above  definitions  for  L,  and  Lq  are  used,   ^T(CJ)^ 
can  be  reduced  to 


+  <Lt^^-  Ci-rH  J  ( I-  P)L  y^ 


(2.32) 


The  operator  Lq  can  be  written  in  terms  of  the  Hamil- 
tonian  for  the  atom  and  electron  subsystems, 

«/       * 

I    =   L,  t    L^    ,  (2.33) 
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vhere  ^ 

(2.34) 
and 

<>  (2.35) 

with  the  assumption  that  the  density  matrix,  ^    ,   Is 
diagonal  In  the  electron  states.  It  follovs  from  symmetry 
that  <L^  >  and  <Lq'^   are  zero.   Hence,  equation  (2.32) 
becomes 

where 


The  matrix  elements  of  ^L^  y   are  specified  In  Section 
(II. 6). 

Five  Identities  derived  In  Appendix  B  that  are 
required  to  further  reduce  JV  (^  )  are  given  belovi 


<^L(»-P)/^>  =  <Z.^(»-P;/^> 


(2.37) 
(2.38) 
(2.39) 
(2.40) 
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i^ir^.(i-p)L,:i"'=  0-p)C*^-<-O  (2.41b) 

vhere  A  is  an  arbitrary  operator. 

Using  the  above  identities,  ^  (6i>)  can  be  reduced  to 

9  12 
an  expansion  in  povers  of  the  operator  L, ,  * 


^«..>-<L,£(^^.)> 


(2.42) 


Introducing  the  unperturbed  resolvent  operator,  K.(6i>), 

defined  by 

-  I 


fC  ito)   2  C'^^  -1-6^ 


equation  (2.42)  becomes 


j.^f^^=  <L,  J(m->('-pH,/>)  > 

-He:/  •  yZt'i.i) 

The  first  term  in  the  above  expansion  becomes 

=  <L,fr,(^^L,>^  (2.44) 

In  this  dissertation  the  above  term  is  the  highest  order 
term  used  in  the  ntmerical  results  presented.  However,  the 
order  of  magnitude  and  the  qualitative  behavior  of  higher 
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order  terms  are  examined  in  Section  (II .7) •  It  Is  shown 
that  the  higher  order  terms  may  become  Important  In  the 
line  center. 
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(II. 6)  MATRIX  ELEMENTS  OF  THE  OPERATOR  ^T(C»>)^ 

To  Obtain  matrix  elements  of  the  tetradlc  ^T(«0)S  , 
the  matrix  elements  of  lJ  and  *f^(lM))   are  required.  Let  the 
Indices M  »^  '   represent  atomic  states  and  «<  »  p  represent 
electron  states.  The  matrix  elements  of  /L?  S  given  In 
equation  (2.43)  are  obtained  by  performing  an  electron 
average  of  equation  (B.9) 

-  ^   ^0-  (B.9) 

The  Boltzmann  factor  f  g  ,  defined  In  equation  (B.4)  Is 
assumed  to  be  diagonal  and  Is  normalized  to  1.  Hence, 
equation  (2.45)  becomes 

(2.46) 
The  matrix  elements  of  ^9^(6i?)  ^  can  be  written  as 
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OMVfjWv 


-  21 


^  '■i(2l"'ir'K^^^P^P. 


(2.47) 


Matrix  elements  of  the  operator  Kq(  Co))  appear  in  equation 
(2*47) •  The  operator  Kq(cJ)  contains  L^t  hence,  to  derive 
an  expression  for  the  matrix  elements  of  Kq(6b»  appearing 
in  equation  (2.47),  the  matrix  elements  of  L^  are  required. 
Combining  equations  (2.46)  and  (B.lO)  one  has 

(2.48) 

To  simplify  the  matrix  elements  of  Kq(^)»  Kq(c>')  is 
expanded  in  powers  of  Lq/«/  and  then  matrix  elements  are 
taken  of  the  individual  terms  in  the  expansion.  In  order 
to  carry  out  this  procedure*  the  matrix  elements  of  L^  are 
required.  One  obtains » 
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M 


yl  L 


vhere  the  results  given  In  equation  (2*48)  have  been  used. 
Resumming  the  above  results  leads  to 

where  the  superscript  m  signifies  that  transitions  between 
initial  and  final  states  due  to  the  ion-field  are  included 
in   ft".  The  above  results  are  inserted  into  the  expression 
for  it^^  (CO)   given  in  equation  (2.47).  Using  the  proper- 
ties of  the  delta  f\mctionSf  we  can  perform  two  of  the 
summations*  The  results  are» 


=  S  (<A<,\  V^  I M w" >  i  ^.  i ^.  -  <vf4l  V^  \vK>  i;^..  ^^.  ) 


44» 
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vThen  the  above  expression  is  multiplied  out,  four  terms 

are  obtained.  Many  of  the  summations  occurring  in  these 
terms  can  be  performed  using  the  properties  of  the  Delta 
function*  The  final  result  after  carrying  out  these  sum- 
mations is  ^  -/ «   .  I  V 


a"V'"' 


H* 


.  J  F„<;.-. ^J i.y> <^,,.^„^  <vWCi^->. 


It"' 


(2.51) 


If  the  transitions  betveen  initial  and  final  atomic 
states  due  to  the  ion-field  are  neglected  and  the  no- 
quenching  approximation  made  for  electrons »  additional 
sximmations  in  equation  (2.51)  can  be  performed.  Under  the 
no-quenching  approximation  (see  Section  1.4),  V     i*»is 
zero  unless  m   and  4\  *  have  the  same  principal  quantum 
number.  The  effect  of  neglecting  these  terms  is  small  for 
initial  and  final  states  that  are  veil  separated  in  energyi 
this  is  the  case  for  the  Lyman-alpha  line.  Hovever*  since 
the  line  profile  is  given  in  terms  of  an  integral  over  all 
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Ion  fields*  one  always  encounters  fields  large  enough  to 
cause  transitions  between  initial  and  final  states*  However » 
the  behavior  of  \Lq(^   )  y    occurring  in  T^^)  is  such 
that  the  Integrand  is  small  for  such  large  fields  and, 
except  for  the  far  wings  of  the  line*  can  be  neglected i   i.e., 
for    |A>)  ^  ISA*  For  the  line  profile  computations 
presented  in  this  dissertation,  the  transitions  between 
initial  and  final  states  arising  from  the  ion  fields  will 
be  neglected* 

It  simplifies  the  derivation  of  ^ (^)  if  the  parabolic 
representation  for  the  atomic  states  is  used*  In  the  para- 
bolic representation  the  z  conponent  of  the  atomic  matrix 
element  R  is  diagonal  with  respect  to  the  principal  quantum 
niimber*  Using  this  representation  and  choosing  a  coordinate 

system  for  the  ion-field  integration  such  that  P,     lies  in 

iz 

the  R  direction,  the  operator  L^  reducei  to 

(2.52) 

Vhen  this  expression  is  used  in  the  expansion  for  Kq(6u)  one 
can  show  that 

^  (2*53) 

When  the  above  foirm  for  K  (^)  is  used  in  the  equation  for 
^  ^^  (pj ) ,  the  matrix  elements  redixce  to  an  expression 
similar  to  those  in  equation  (3*119)  in  reference  9i 
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mAM 


'i^ti  \  m"^ 


(2.54) 

nhere  the  nm  superscript  indicates  that  the  transitions 
between  initial  and  final  states  due  to  the  ion-field  are 
not  included.  Note  that  k^  contains  the  ion-field,  vhereas 
the  corresponding  expression  in  reference  9  does  not  contain 
the  ion-field. 

Computational  forms  for  equation  (2.54)  are  derived  in 
Appendix  C.  If  V^^  is  taken  to  be  the  full  Coulomb  inter- 
action and  if  symmetry  effects  in  P  are  neglected,  one 

obtains 

oo      

-tf"   o      *- 

/,    __L_  ^_L_ 

X  V"   l-t-  (\«-/iiTr/  f  (2.55) 
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vhere  ^  "  fc  T  »  ^MiA   ^"  defined  In  Appendix  D,  and  A^ 
is  the  Debye  length.  The  above  expression  is  the  result 
of  evaluating  the  first  term  in  equation  (2 •54).  If  lover 
state  interactions  are  neglected,  the  other  terms  do  not 
enter  into  the  computations.  In  the  dipole  approximation 
the  lover  state  interaction  is  identically  zero  for  the 
Lyman-alpha  transition.  However,  vhen  the  full  Coulomb 
interaction  is  used,  the  lover  state  interaction  is  not 
zero.  Computations  for  the  Coulomb  case  were  performed  in 
which  the  lover  state  interaction  vas  Included  i  the  results 
show  that  these  lower  state  interactions  produce  less  than 
a  5  percent  effect  on  the  line  profile  in  the  temperature 
density  cases  examined.  However r  these  interactions  may 
have  a  larger  effect  for  the  higher  series  transitions. 

From  the  above  analytic  form  it  appears  that  vhen  the 
integration  variable  approaches  zero,  convergence  problems 
arise  in  the  line  center.  This  votO-d  in  fact  be  the  case 
if  electron  correlations  were  neglected  in  the  electron 
distribution  function.  However,  the  second  term  in  brackets 
contains  electron  correlations  and  provides  convergence  in 
the  line  center.  When  electron  correlations  are  neglected, 

the  resulting  divergence  in  ?^  (£•> )  produces  a  dip  in  the 

9  35 
center  of  the  computed  profile. 
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(n.7)  HIGHER  ORDER  TERMS  IN  5^(^) 

Most  of  the  Interactions  betveen  the  radiator  and 
perturbing  electrons  are  veak  because  the  probability  is 
small  that  electrons  get  close  enough  to  the  radiator  to 
cause  a  strong  collision  (see  Section  II •2).  Silnce  ^  {<t)) 
is  expanded  in  terms  of  the  coupling  constant,  it  is 
expected  that  the  perturbation  expansion  for  ^{<a))   vill 
converge  rapidly  for  veak  interactions  or  veak  coupling. 
Hotrever»  in  vhat  follovs*  it  is  shoyn  that  higher  order  terms 
may  become  appreciable  for  small  frequencies i  the  pertur- 
bation expansion  breaks  dovn  as  A<b> — ^0  (:t-*'cao). 

The  computation  of  higher  order  terms  Is  not  tractable. 
Hence,  only  an  approximate  analysis  of  the  behavior  of  a 
higher  order  term  in  the  expansion  for  ^(Ca))  is  considered. 
The  procedure  folloved  is  a  generalization  of  the  method 
used  in  obtaining  an  expression  for  the  first  term  (see 
Appendix  C). 

The  n  K  2  terms  in  the  expansion  vanish  because  of  the 
symmetry  of  the  electron  distribution  function  and  the 
radiator-electron  interaction.   Thus,  the  next  contributing 
term  is  the  n  ■  3  term.  There  are  (n  •>>  1)  «  16  such  terms. 
A  typical  term  has  the  form 
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(2.56) 

The  other  fifteen  terms  are  obtained  by  permutations  of  the 
indices  in  the  above  term* 

For  the  purpose  of  estimating  the  magnitude  of  the  , 
above  expression  the  no-quenching  approximation  is  made  and 
the  ion-field  splitting  of  the  atomic  levels  is  neglected, 
under  these  conditions  the  operator  k.  can  be  written  as 

Ob 

^t^P       V  (2.57) 

Inserting  these  results  into  equation  (2.56)  yields 
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The  terms  ^ju<'  *£'«  ^^^   change  in  energy  of  the 
perturbing  electron  during  a  collision.  At  the  line  center 
^•(*  Is  snail  because  changes  In  the  energy  of  the  atomic 
electron  are  small  and  there  are  no  other  significant 
mechanisms  in  this  theoretical  treatment  vhich  could  cause 
a  large  change  in  the  perturbing  electrons*  energy*  Also, 
the  terms  Co)  «  ^^  are  small  in  the  line  center.  Thus* 
at  the  line  center  the  exponentials  do  not  oscillate  rapidly 
as  they  do  in  the  vlngs  of  the  line*  Hence,  it  is  expected 
that  higher  order  terms  may  become  large  in  the  line  center 
but  not  in  the  line  vlngs  since  the  oscillating  nature  of 
the  exponentials  assures  that  the  above  term  is  small  in 
the  line  wings. 

If  the  procedure  used  in  Appendix  C  is  generalized,  an 
Integral  expression  for  the  n  «  3  term  in  ^  (<^ )  can  be 
obtained.  The  following  multi-dimensional  Integral  is 
encountered I 

A </\r'><  f- 1*><* I  f'?i,?"[ury^ur- 1 r"><r"'  1  ^y 

Tfhere  the  meaning  of  the  various  terms  in  this  equation  and 
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in  subsequent  equations  in  this  section  are  explained  in 
Appendix  C.  After  transforming  to  relative  coordinates  and 
performing  all  possible  Integrals  analytically,  one  obtains 

'<fw^^^iv"  ^'^'-*'^W  ('^-?'Oi..v  (^"> 

(2.59) 

The  electron  correlations  provide  a  convergence  factor 
for  the  q  integration  as  they  do  in  the  first  term,  n  ■  1, 
of  the  expansion  (see  Section  IT .6).  However,  the  q*  and 
q**  integrations  do  not  have  suitable  convergence  factors  for 
small  frequency  separations*  The  results  of  performing  the 
q**  Integration  vill  be  similar  to  the  results  obtained  in 
the  n  a>  1  term  for  1  ■  j  (see  equation  C.23,  Appendix  C). 
This  factor  diverges  in  the  line  center.  Thus,  the  q*  and 
q**  integrations  yield  results  which  diverge  in  the  line 
center.  Hence,  the  higher  order  terms  are  appreciable  in 
the  line  center. 

The  above  arguments  hold  whether  the  full  Coulondb 
interaction  or  the  dipole  approximation  is  used.  However, 
there  are  differences  between  the  two  treatments  which  are 
considered  next. 

In  the  Lyman-alpha  transition,  four  matrix  elements 
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appearing  in  the  first  order  term  (n  »  1 )  for  ^  (^ )  are 
found  to  account  for  most  of  the  broadening*  When  these 
elements  are  calculated  using  the  full  Coulomb  interaction 
they  are  larger  than  the  corresponding  terms  in  the  dipole 
approximation*  The  ratio  of  the  dipole  approximation  to 
the  full  Coulomb  interaction  treatment  for  one  of  the  most 
sigaificant  matrix  elements  is  shown  in  Figure  2*  The  ratio 
is  alvays  greater  than  one  for  the  various  strong  collision 
cut-offs  considered*  Therefore*  the  perturbation  expansion 
for  ^  ((a>)   using  the  full  Cotilorab  interaction  vill  break 
dovn  somevhat  sooner  as  ^^  approaches  zero  than  it  will 
when  the  dipole  approximation  is  employed*  This  behavior 
may  partially  explain  the  differences  between  the  two 
treatments  in  the  line  center* 
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CHAPTER  III 
QUALITATIVE  ANALYSIS  OF  TOE  LINE  SHAPE  EXPRESSION 

(III.I)  INTRODUCTION 

The  mathematical  aspects  of  the  line  shape  problem 
necessary  to  Interpret  the  theoretical  results  are  analyzed 
In  this  chapter*  Expressions  for  the  line  shape  are  derived 
In  vhlch  the  effects  of  various  approximations  can  be  easily 
related  to  the  computed  profiles.   In  certain  cases, 
functional  aspects  of  the  mathematical  expressions  are 
examined  rather  than  the  complete  expressions  used  In  the 
computations. 
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(III. 2)  THE  OPERATOR  O^-  <T(C«>)> 
To  compute  a  line  profile  the  expression 

must  be  evaluated.  This  Implies  that  the  matrix  elements 
of  the  operator  Ci>  -  <T(a))>  must  be  calculated.  Due  to 
the  Llouvllle  operators  occurring  In  <T(^ )^  ,  these  matrix 
elements  will  be  high  order  tensors  requiring.  In  general, 
four  or  more  Indices.  The  matrix  elements  of  the  Llouvllle 
operator,  L^,  and  of  the  effective  Interaction  operator,  ^  (<k) ), 
are  given  In  equations  (2.48)  and  (2.55)  respectively.  When 
these  are  used  In  the  expression  for  (<ti  -  ^T(c«i)^  ),  one 
obtains  the  following i 

(3.1) 
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vhere  C.  and  C2  are  constants.  As  indicated  in  Appendix  C, 
the  function  v/  (^)  depends  on  the  ion-field  shifted  frequency. 

The  first  bracketed  term  in  the  above  expression  repre- 
sents the  frequency  separation  from  the  line  center. 
Mathematically  it  can  be  interpreted  as  the  independent  vari- 
able. 


The  second  term  contains  the  broadening  effects  of  the 
static  ions.   f{^^,    is  restricted  to  lower  states  and  P     , 
is  restricted  to  the  upper  states.  Hence»  lover  state 
interactions  vith  the  ion-fields  are  included.  The  lower 


state  ion-field  broadening  is  identically  zero  for  the  Lyman- 
alpha  line  in  hydrogen.  However,  it  may  become  important 
in  the  Balmer  beta  transition  having  as  much  as  a  10  percent 
effect  in  certain  cases. ^^ 

The  last  two  terms  in  equation  (3.1)  represent  the  effects 
of  electron  broadening.  In  the  relaxation  theory  these  terms 
are  not  symmetric  in  ^  -  <a>  ..  and  are  solely  responsible  for 
a  red  asymmetry  in  the  line  center  that  changes  to  blue  in 
the  far  wings  of  the  line.  The  red  to  blue  cross  over  is 

quite  apparent  for  the  Lyman-alpha  transition  in  ionized 
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helium.   '    It  also  appears  in  the  Lyman-alpha  profile  for 

hydrogen  but  only  in  the  far  wings  of  the  line,  ^^  y     8A. 

As  discussed  in  Section  (II. 6),  lower  state  perturbations 
due  to  electron  interactions  are  contained  in  the  last  term 
in  brackets  in  equation  (3.1). 

In  reference  9  it  is  shown  that  for  the  Lyman-alpha 
transition  only  the  V"  ■  \r*  "  1  terms  contribute  to  the  line 
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profile.  In  the  no-quenching  approximation  the  states  M.M 
are  restricted  to  initial  states.  Thus,  for  the  Lyman-alpha 
transition  only  two  Indices  are  required  to  specify  the 
elements  in  the  tetradic  to-  ^T(Ci>)^  •  That  is  (  ^  - 
^T(u»)^  )  has  a  convenient  matrix  representation.  Further- 
more, since  a  hydrogenic  state  having  principal  quantum 
number  n,  has  n  degenerate  levels,  the  matrix  representation 
vill  be  x?   X  rrx   it  will  also  have  a  block  diagonal  form  as 
will  be  8hot.m  in  the  next  section. 
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(III. 3)  THE  PARABOLIC  REPRESENTATION 

The  Schroedlnger  equation  vlthout  an  external  field  is 

separable  in  the  parabolic  coordinate  system.  The  resulting 

38 
vave  functions  are  given  by 

//  l^*U(  1 

^M^M^M  U     J     /-     >H 

(3.2) 

vhere  ^^^^  aire  the  parabolic  quantum  numbers  associated 
with  the  atomic  state >  4:M\  ,     The  other  symbols  are  defined 
in  Appendix  D.  It  should  be  emphasized  that  the  parabolic 
wave  functions  given  above  are  not  '*Stark**  wave  functions 
since  they  are  not  solutions  to  the  Schroedinger  equation 
for  a  system  consisting  of  a  hydrogen  atom  in  an  external 
field. 

The  parabolic  representation  is  used  for  mathematical 
convenience.  For  a  given  principal  quantum  number,  the  t 
component  of  the  dipole  matrix  element  is  diagonal  in  this 
representation.  Hence,  the  atomic  density  matrix  occurring 
in  the  line  profile  expression  is  diagonal}  the  second  term 
in  brackets  in  equation  (3.1)  is  also  diagonal,  and  many  of 
the  other  matrices  reduce  to  a  block  diagonal  form.  These 


49 
slmpllflcatlona  are  particularly  convenient  for  computational 

purposes* 

In  the  full  Coulomb  Interaction  case*   the  term     C^^ii  Xu"^' 

occurring  In    £i> '"  ^T  (<a? )  ^      contains  a  delta  function  factor 

of  the  form 


-    /„ 


-m'^m!'  '%'"*!'  ^'^^r^* 


vhere  m  is  the  magnetic  quantum  number.  Thus*  ^mm"  ^M"4t' 
is  diagonal  in  m  when  the  full  Coulomb  Interaction  is  usedo 

In  the  dlpole  approximation  the  third  term  in  the 
expression  for  ^  -^tC^)"^,  given  in  equation  (3*1)^  becomes 


with  the  results  of  Appendix  F,  the  dot  product  (\uu'^  *    m'*4A* 
is  found  to  contain  delta  function  products  of  the  form 


MM         "^M  ^ 


^m' 


Thus  in  both  the  dlpole  interaction  and  the  full  Coulomb 
interaction*  the  matrix  (  (^   -  ^(^)^)    has  the  form 
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(<y-  <T/<t>)>^ 
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where  the  parabolic  states  are  labeled  as  followst 


ll>   a 

|100> 

12)  B 

|20-1> 

|3>  S 

1200) 

147  S 

1210  > 

157  3 

1201  >  . 

Blocks  vhich  arise  frcm   the  same  magnitude  of  m  but 
different  signs  are  mathematically  equivalent  and  can  be 
shovn  to  give  identical  contributions  to  the  line  profiles 
thus*  the  22  block  contributes  the  same  as  the  55  block* 
This  follows  from  the  form  of  the  wave  function  given  in 
equation  (3.2)  t  the  wave  function  depends  on  the  sign  of  m 
only  through  the  term  exp(im^)  ^ich  does  not  affect  the 
magnitude  or  sign  of  any  of  the  matrix  elements  (see 
Appendices  D  and  F)  • 

Following  a  similar  procedure  for  the  Balmer  beta 
transition,  it  is  found  that  the  orders  of  the  matrices 
encountered  range  from  1  to  4« 
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(III. 4)  THE  LYMAN-ALPHA  TRANSITION 

From  the  previous  section.  It  Is  clear  that  the  symmetry 
Of  the  atomic  states  In  the  parabolic  representation  greatly 
simplifies  the  computational  effort.  Hovever.  there  still 
remains  the  problem  of  reducing  the  line  shape  expression 
given  In  equation  (2.25)  to  quadratures.  This  section 
outlines  the  procedure  for  reducing  the  Lyman-alpha  line 
shape  expression  to  a  convenient  calculable  form. 

The  matrix  elements  of  the  operator  (C*)  -  <t(U))>  )  are 

^f<    S    ^  ^  /  ^'*^'^*  fswA^y!"^  q ,  / 

8      "       ^ 


(3.3) 


vhere 


€      »  Ion  mlcrof leld 


W 


MM 


electron  collision  term. 


In  the  dlpole  approximation  \A^  i  Is  given  by 

«here   Z^-,   contains  the  lon-f leld.  Using  the  results 
derived  In  Appendix  F  «e  find  that  the  non  zero  matrix 
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elements  of   NA/  i  are 


XI /A 


4- J       i3  4- 


(3.5) 


T-rhere   A^O^   =  A^^  =  4<iJ.,  is  the  unshlfted  frequency 
separation.  The  corresponding  elements  of  (  6u  -  4^T(CLi)y)   are 


31-   - 


(3.6) 
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vhere  2-33  -  38^.  2r^^  m   -sa^,  and  the  definition  of  c 
follows  from  equation  (3.3), 


All  other  matrix  elements  are  zero  or  do  not  enter  Into  the 
computations  because  of  the  restriction  that  X  have  matrix 
elements  only  between  the  Initial  and  final  states. 

The  resolvent  operator,  K(<a>),  appearing  In  the  line 
shape  expression  Is  the  Inverse  of  ^  -  <;t(6o)>  ,  por  block 
diagonal  matrices  one  has 


J_Aft-5f\   ' 


vhere  D  Is  the  determinant  of  the  2  x  2  matrix. 

The  determinant  of  the  2  x  2  block  reduces  to 

t^hlch  can  be  rewritten  In  the  following  convenient  formi 


(3.7) 
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■where  A  and  B  are  complex  coefficients  that  depend  on  the 

matrix  elements  of  W*  A  and  B  are  small  compared  to 
J\  AdJ^/  -hr   except  in  the  line  center  (see  Chapter  V  for 
actual  nimierical  comparisons)*  Hence »  the  determinant 
provides  a  resonant  denominator  iihich  has  a  minimum  at  ion- 
fields  given  by 

(3.8) 

The  function  J((t;  ,  g  )  used  to  compute  the  line  profile  (see 
equation  2.15)  will  be  peaXed  at  values  of  the  field  given 
in  the  above  equation.  J(Ci;  ,€)   can  be  interpreted  as  a 
line  profile  from  a  radiating  atom  that  is  placed  in  a 
constant  ion-field  and  perturbed  by  the  electrons*  This 
interpretation  is  useful  since  considerable  physical  insight 
into  the  effects  of  electron  broadening  can  be  gained  by 
studying  the  function  J(^,€)   rather  than  the  full  line 
profile* 

For  transitions  more  complicated  than  the  Lyman-alpha 
profile*  there  win  be  additional  resonant  denominators 
corresponding  to  the  various  Stark  shifted  levels*  For 
example,  in  the  Balmer  beta  line  there  are  three  different 
z  matrix  elements  resulting  in  three  distinct  resonant 
denominators* 

To  compute  a  line  shape,  the  imaginary  part  of  the 
following  trace  is  required 
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The  matrix  A    occurring  in  the  above  expression  is  restricted 
to  have  elements  only  betveen  initial  and  final  states.  The 
matrix  elements  of  the  x,  y,   and  z  components  are  computed 
in  Appendix  F« 

In  the  no-quenching  approximation »  the  atomic  density 
matrix  f   can  be  vritten  as, 


e 

Using  the  above  results  the  trace  over  atomic  states 
can  be  performed.  The  results  arei 

(3.9) 

The  leading  term  in  brackets  does  not  have  a  strong  ion- 
field  dependence  in  the  dipole  approximation  and  is  field 
independent  in  the  full  Coulomb  interaction  treatment.  In 
the  dipole  approximation,  this  term  can  be  interpreted  as 
representing  the  effects  of  electron  broadening  on  the  initial 
levels,  I  3^  and   |4^  .  In  the  dipole  approximation  the 
K^2   ^®i^  depends  on  the  dipole  matrix  elements  having 
different  magnetic  (m)  quantum  number  as  can  be  seen  by 
examining  the  vVjj  matrix  element  given  in  equation  (3.4). 
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Hovever,  the  full  Coulomb  interaction  does  not  mix  states 
having  different  m  values i  as  a  consequence,  the  leading 
term  is  field  independent* 

In  both  the  dlpole  approximation  and  the  full  Coulomb 
interaction  case*  the  K..  contribution  is  significant  only 

in  the  line  center*  For  a  temperature  of  20,000  K  and  an 

17   3 
electron  density  of  3  x  10  /cm  the  K22  contribution  for 

o 

the  dlpole  approximation  is  slgnlflcsmt  from  0-1  A  from 
the  line  center*  Whereas,  for  the  full  Coulomb  Interaction 

o 

treatment  the  range  of  Importance  is  from  0  >  *8  A*  The 
net  effect  of  this  slightly  different  behavior  is  to  produce 
a  shoulder  in  the  line  shape  vhen  the  full  Coulomb  lnter> 
action  is  used*  It  is  interesting  to  note  that  a  similar 
behavior  is  encountered  in  the  dlpole  approximation  at 
higher  temperatures  in  hydrogen  and  in  the  Lyman -alpha 
transitions  in  ionized  helium* 

By  taking  the  VsZ-f-O  limit  in  equation  (2*54)  a  static 
wing  formula  can  be  obtained*  This  is  the  subject  of  the 
next  chapter* 


CHAPTER  IV 
ASYMPTOTIC  WING  FORMULA 

(IV. 1)  INTRODUCTION 

The  purpose  of  the  present  chapter  Is  to  derive  an 
expression  that  accurately  describes  the  wings  of  the  profile. 
In  Section  (IV. 2)  a  heuristic  derivation  of  a  static  vlng 
formula  is  presented.  The  restating  formula  is  shown  to 
be  the  leading  term  in  a  rigorous  perturbation  expansion 
presented  in  Section  (IV. 3).  Higher  order  terms  in  this 
expansion  are  shown  to  Include  the  effects  of  electron 
dynamics • 


57 


58 


(IV. 2)  HEURISTIC  DERIVATION  OF  A  STATIC  WING  FORMULA 


In  the  derivation  considered  In  this  section,  it  is 
assumed  that  the  only  broadening  mechanisms  are  static  fields. 
Under  this  assumption  a  wing  formula  can  be  obtained  from 
the  line  profile  expression  (equation  3.9)  presented  in 
Chapter  III  by  the  follovlng  procedure  i  1.  Replace  all 
static  ion-fields,  £^.  ,  by  the  total  static  field, 
^j  (     ^y"  ^i*  ■*■  ^  ^'   ^*  Take  the  limit  as  the  electron 
dynamic  terms  go  to  zero*  In  this  limit  the  line  profile 
becomes 


+      ^/OP«0 


(4.1) 


vhere  P,j,(  ^   )  is  the  microfield  probability  function  for 
the  combined  electron- and  ion-fields,  defined  as  £_C]^(  £'<p  )  • 
Since    P(^)w,        f^(-  ^  )  (see  Section  III. 4)  and 
P—  (  ^•.  )  is  symmetric 
equivalent  to 


14 


in 


£   ,  the  above  expression  is 


STATIC 


=   A/^,C^.)Pt^T) 


+  ^K 


4a 


^'= 


(4.2) 
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The  atomic  density  matrix  element  (^3  evaluated  at 

^T  ~  3tf  <  €    ^^   equal  to  1/e  where  e  Is  the  natural  loga- 
rithm base*  Thus 

(4.3) 

Determining  the  mlcrofleld  distribution  function  for 
the  total  field  from  first  principals  Is  non  trivial.  As  a 
first  approximation  a  convolution  of  the  electron  mlcrofleld 
vlth  an  Ion  mlcrofleld  Is  used.  The  procedure  for  obtaining 
the  total  mlcrofleld  In  terms  of  a  convolution  follows. 

The  mlcrofleld  Is  defined  as 


QUr)    =  \.(^.  ^iK-^i-t^) 


(4.4) 


By  writing  the  density  matrix  In  factored  form,  one  can  show 

39 
that  the  above  reduces  to 

=  J  A'Ti,  ( p.  i,  (t-  £,  )]\  (f^  S  ( vt-t,-)) 

J  (4.3) 

where  the  properties  of  the  Dlrac  delta  distribution  have 
been  used  to  obtain  the  last  expression  which  Is  a  convolution 
of  mlcroflelds. 
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It  should  be  enphaslzed  that  the  above  expression  for 
the  total  mlcrofield  Is  an  approximation t   in  order  to 
obtain  the  factored  form  for  the  density-  matrix,  the 
electron-ion  interaction  vas  taken  into  account  by  an 
effective  potential  vhich  depends  only  on  the  ion  coordi- 
nates (see  Section  II. 1  and  Section  1.4) •  Hence,  the 
electron-ion  interaction  is  contained  only  in  the  ion 
microfield.  It  is  expected  that  this  treatment  of  the 
electron-ion  interaction  should  be  good  since  the  mobility 
of  the  electrons  is  much  larger  than  the  ion  mobility 
resulting  in  a  shielding  of  the  ion-ion  interaction.  To 
examine  the  effects  of  shielding  on  the  electrons-electron 
interaction,  line  profiles  obtained  from  convolutions  of 
various  combinations  of  lov-and  high-frequency  microfields 
are  compared  in  Chapter  V« 
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(IV. 3)  THEORETICAL  DERIVATION  OF  A  WING  FORMULA 

In  the  line  wings  the  time  of  Interest  Is  small t  hence, 
a  perturbation  expansion  In  which  the  time  Is  an  expansion 
parameter  may  prove  to  be  useful  In  the  line  wings.  Such 
an  expansion  Is  derived  In  this  section  and  carried  to 
second  order.  The  procedure  for  generating  higher  order 
terms  follows  from  a  straightforward  extension  of  the 
method  used  to  calculate  the  second  order  term. 

The  Hamlltonlan  appearing  In  the  function,  <$(  X),   can 
be  written  as 


where 


H  -  H,  +  H^ 


(4.ea) 


W.  =  W,+  Hi    +V^,-. 


(4.6b) 


Using  the  above  definitions  the  time  development 
operator  T(X')  appearing  in  25  (JT  )  can  be  written  as 

where  exponents  containing  third  and  higher  order  terms  have 


62 
been  neglected i  i.e.* 


.e,y^(-x'"C/<,;-W,H,>.])  ^   1.  forn>2. 
Similarly 

Substituting  the  above  results  in  to  ^(^)  yields 

(4.8) 

It  is  convenient  to  regroup  the  above  equation  into 
the  following  formi 


(4,9) 
vhere  again  terms  such  na  x/Lft*(C^:t  )  have  been  set  equal  to  one. 
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The  following  operator  expansion   is  used  to  further 
reduce  $  (  ^)  S 

•  (4.10) 

vhere  A  and  B  are  operators  defined  as 

(4.11) 
and  then  applying  the  above  to  equation  (4.9)  for  ^(jC)yL«ld8 

+   ^^*'). 

(4.12) 

Clearly  the  above  expression  vill  break  down  for  large 
times.  Hence*  as  smaller  frequency  separations  are  approached, 
higher  order  terms  vill  become  increasingly  important.  It 
is  expected  that  electron  dynamics  will  certainly  be 
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significant  for  frequencies  less  than  the  electron  plasma 

frequency!  thus*  for  the  region  of  the  profile  corresponding 

to^4)<iA^  terms  linear  in  :C  and  perhaps  even  higher  order 

terms  ^rill  be  required  to  adequately  represent  the  intensity. 

Based  on  comparisons  vith  experimental  data  presented  in 

Chapter  V»  this  behavior  indeed  appears  to  be  verified. 


CHAPTER  V 
ANALYSIS  OF  THE  THEORETICAL  RESULTS 

(V.l)  INTRODUCTION 

An  analysis  of  the  theoretical  computations  is 
presented  in  this  chapter.  The  theoretical  results  are 
compared  with  the  experimental  data  and  with  other 
theoretical  computations. 
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(V.2)  THE  FUNCTION  J(U>,S  ) 

The  line  profile  involves  an  integration  over  £  of 
the  product  of  P(  6  )  times  J(iOf  £  ).  The  region  of  the 
profile  in  ^ich  various  field  dependent  effects  are 
important  can  be  determined  by  an  analysis  of  J(£^  «9  )• 

Ji(jit€  )   consists  of  a  field  independent  term  plus  a 
field  dependent  termi 

J(o*  »€  )  ■  J(«*)  +   ^(W»e).  (5.1) 

In  the  case  of  Lyman-alpha  radiation^  a  graph  of  fU  (co,^) 
vs*  £*  vill  have  one  sharp  peak  due  to  the  resonance 
denominator  discussed  in  Section  III .4  (see  Figure  3). 

A  slightly  more  complicated  result  is  found  vhen 
^  (<A> ,  £  )  vs.  E    is  plotted  for  the  Lyman-beta  transition  in 
hydrogen.   In  this  case,  the  ion-field  removes  the  degeneracy 
in  the  upper  level*  splitting  it  into  three  distinct  levels* 
tvo  of  vhich  are  shifted.  Hence,  the  graph  of  the  field 
dependent  part  of  J(Cu  ,£  )  vill  have  tvo  peaks.  /^  {<** »  i  ) 
for  the  Lyman-beta  transition  is  plotted  in  Figure  4.  For 
£   a  0  both  peaks  coincide  auid  occur  at  £  a  o.  As  the 
frequency  increases  tvo  peaks  appear. 

A  fev  general  characteristics  of  A  (oj,£  )   and  some  of 
the  resulting  implications  are  listed  belov. 


Figure  3.  Field  dependent  part  of  the  fvinction  J(W,  £  ) 
for  the  Lyman-alpha  profile. 
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Figure  4.  Field  dependent  part  of  the  function  J(u>  ,£  ) 
for  the  Lyman-beta  profile. 
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1)  For  large  frequency  separations  from  the  line 

center f  most  of  the  contribution  to  the  field  integration 
appearing  in  the  line  shape  expression  arises  from  large 
values  of  the  field* 

2)  The  relative  maximum  of  JD(u>>£)  increases  as  the 
frequency  increases  i«hile  the  half -width  decreases.  This 
suggests  that  replacing  the  peaks  in  the  function  J|(cAj»g  ) 
by  delta  function  distributions  should  be  accurate  in  the 
line  wings*  This  limiting  behavior  was  used  in  the  heuristic 
derivation  of  the  static  wing  formula  discussed  in  Section 
(V.2\ 

3)  The  maximum  value  of  J  (^u*  £  )  decreases  rapidly  as 
the  frequency  increases*  but  the  rate  of  decrease  diminishes 
for  larger  values  of  the  frequency. 

4)  For  every  value  of  ^CO   >  ^(<«>»6.)  will  approach  a 
fairly  constant  value  for  field  values  greater  than  those 
for  which  ft  (^»  £  )  peaks. 

From  the  behavior  of  J(Co  >  g  )  it  is  easy  to  determine 
the  range  in  frequency  for  which  various  field  dependent 
effects  are  important.  For  example,  a  field  dependent  effect 
that  is  directly  proportional  to  the  strength  of  the  field 
will  be  most  important  in  the  line  wings.  One  such  effect 
which  exhibits  this  dependency  arises  from  the  field  depen- 
dence in  the  atomic  Boltzmann  factors.  This  is  the  subject 
of  the  next  section* 
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(V.3)  THE  ION-FIELD  DEPENDENT  ATOMIC  DENSITY  MATRIX 

The  Stark  shift  In  energy  levels  caused  by  the  ion-field 
is  frequently  neglected  in  the  atomic  density  matrix  occur- 
ring in  the  line  shape  expression.  Neglecting  the  field 
assumes  that  the  density  of  states  is  equal  for  all  atomic 
levels  corresponding  to  a  given  principal  quantum  number. 
The  ion-field  removes  the  degeneracy  of  the  excited  states 
via  the  Stark  effect  yielding  2n  -  1  distinct  energy  levels 
or  Stark  levels  tfhere  n  is  the  principal  quantum  number* 
The  field  dependent  atomic  density  matrix  -ueights  Stark  states 
that  are  lover  in  energy  more  heavily  than  states  of  higher 
energy.  Since  transitions  from  lover  energy  initial  states 
give  rise  to  the  red  ving  of  the  line  profile,  inclusion  of 
the  ion-field  in  the  atomic  density  matrix  vill  cause  a  net 
red  asymmetry.  To  determine  vhether  this  effect  is  uniform 
over  the  entire  profile  one  can  analyze  the  fxmction  ^  (u)  ,€  ), 
In  the  previous  section,  curves  that  demonstrate  the  behavior 
of  ^(t*^,  ^)   vere  presented.  For  the  Lyman-alpha  transition 
the  peak  in  ^(£^»^  )  occurs  at  a  field  value  given  approxi- 
mately by 

Thus,  for  large  frequency  separations,  J(<«>,^)  peaks  at 
large  values  of  the  ion  microfield.  The  exponent  appearing 
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In  the  Boltzmann  factors  In  the  atomic  density  matrix  is 

directly  proportional  to  the  lon-fleld  and  hence,  the  larger 

the  field 9  the  greater  the  effect  in  the  density  matrix* 

Thust  the  effect  of  Including  the  field  In  the  atomic 

density  matrix  Is  relatively  more  Importzmt  In  the  wings  of 

the  line  than  It  Is  In  the  line  center*  Profiles  calculated 

vlth  and  vlthout  an  lon-fleld  dependent  atomic  density  matrix 

are  given  In  Figures  5  and  6  respectively.  From  these  curves 

It  is  demonstrated  that  Including  the  Ion-field  can  produce 

0 

a  noticeable  effect  for  wavelengths  greater  than  4  A. 

The  dependence  of  the  atomic  density  matrix  on  Lempers- 
ture  and  on  the  magnitude  of  the  dlpole  matrix  element  should 
also  be  considered*  The  lon-fleld  dependent  factor  appearing 
In  the  exponent  of  the  atomic  density  matrix  Is  proportional 
to  g  R^  /-A J"  ^ere  T  Is  the  temperature  and  R^  Is  the  z 
component  of  the  atomic  dlpole  operator*  Thus,  at  low 
temperatures,  the  Importance  of  the  Boltzmann  factor 
Increases*  Also,  for  higher  series  transitions  R^  is  larger 
and  the  effect  may  become  large*  Since  radiation  from  higher 
series  transitions  will  be  of  longer  wavelength,  the  corres- 
ponding frequency  separations  will  be  smaller*  For  small 
values  of  the  frequency  separation  the  field  Integration 
over  ^{(ti,^)   contributes  significantly  only  for  small  values 
of  the  field*  Thus,  for  higher  series  transitions,  the 
product  £'R^  is  not  as  large  as  one  might  otherwise  expect* 

For  the  Lyman-alpha  transition  In  Ionized  helium,  the  z 
component  of  the  atomic  dlpole  matrix  element  Is  smaller  by  a 


Figure  5.  Lyman-alpha  profile  using  a  field  dependent 
atomic  density  matrix* 
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Figure  6,  Lynam-alpha  profile  neglecting  the  field 
dependence  In  the  atomic  density  matrix. 
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factor  of  t-wo  compared  to  the  matrix  elements  in  hydrogen* 
Al80»  the  helium  profiles  presented  in  this  dissertation 
correspond  to  temperatures  that  are  about  tvice  as  large  as 
those  considered  in  the  hydrogen  line  calculations.  ThuSf 
the  exponent  appearing  in  the  atomic  density  matrix  is 
decreased  by  a  factor  of  approximately  4,   and  the  effect 
of  the  ion-field  in  the  density  matrix  can  be  neglected. 
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(V.4)  THE  FULL  COULOMB  INTERACTION 

One  of  the  approximations  used  in  most  treatments  of 
the  line  broadening  problem  is  the  assumption  that  the 
interaction  between  the  perturbing  electrons  and  the 
radiating  atoms  is  a  dipole  interaction.  Such  an  approxi- 
mation provides  an  accurate  description  of  the  interaction 
for  large  distances  between  the  atom  and  perturbing  electron* 
To  achieve  10  percent  agreement  between  the  exact  interaction 
and  the  dipole  approximation  for  all  orientations  of  an 
electron-proton  dipole,  the  distance  between  the  dipole  and 
the  point  at  which  the  potential  is  measured  must  be  greater 
than  10  times  the  length  of  the  dipole.   For  excited  states 
resulting  in  Lyman-alpha  radiation,  the  average  distance 
between  the  nucleus  and  atomic  electron  is  4  Bohr  radii. 
Hence*  to  achieve  at  least  10  percent  agreement  between  the 
dipole  approximation  and  exact  interaction,  the  distance 
between  the  perturbing  electron  and  radiating  atom  should  be 

greater  than  40  Bohr  radii.  The  strong  collision  cut-offs 

1  ?  9  16  21 
used  in  various  theories  '  »  »   »   for  the  Lyman-alpha 

transition  are  less  than  15  Bohr  radii  which  is  in  a  region 

vhere  the  dipole  approximation  is  not  highly  accurate. 

These  conditions  provided  the  motivation  for  developing  a 

line  broadening  theory  which  uses  the  full  Coulomb  potential 

to  describe  the  atom-perturbing  electron  interaction. 
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Hovever,  it  should  be  emphasized  that  even  in  the  full 
Coulomb  treatment,  a  strong  collision  cut-off  should  be 
used  to  approximate  symmetry  effects  usually  neglected  in 
the  electron  density  matrix  (see  Section  II •2). 

One  effect  of  including  the  full  Coulomb  interaction 
in  this  second  order  theory  is  to  produce  a  shoulder  at 
about  0.5  A  in  the  Lyman-alpha  profile.  A  shoiader  is  also 
obtained  in  the  dipole  approximation.  Although,  it  is  not 
as  pronounced  as  in  the  full  Coulomb  treatment.  In  the 
parabolic  representation  the  dipole  matrix  elements  for  the 
initial  states  have  non  zero  x  and  y  components  between 
states  having  different  magnetic  quantum  numbers  (see 
Appendix  F) .  Hovever,  matrix  elements  of  the  full  Coulomb 
interaction  are  diagonal  with  respect  to  magnetic  quantum 
number  (see  Appendix  D) .  Also,  the  magnitude  of  the  matrix 
elements  of  W  are  different  for  the  two  approximations.  The 
matrix  elements  which  appear  in  the  function  jiig))   in  the 
full  Coulomb  interaction  treatment  are  smaller  than  the 
corresponding  elements  in  the  dipole  approximation  and 
are  primarily  responsible  for  the  more  pronounced  shoulder. 
The  numerical  values  for  the  W  matrix  for  the  two  cases  are 
compared  in  Figure  7  for  various  strong  collision  cut-offs. 
The  computed  line  shape  results  for  the  full  Coulomb  treat- 
ment are  presented  in  Figure  8. 

With  either  the  full  Coulomb  interaction  or  dipole 
approximation,  the  line  profile  can  be  written  as 

I(a>)  «  j.(oJ)  +   ^p(e)  ^(^.£)^6.       (5*3) 
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Figure  8«  Lyman-alpha  profile  using  the  full  Coulomb 

interaction.  Solid  line  does  not  use  a  strong 
collision  cut-off.  Dashed  line  uses  a  strong 
collision  cut-off. 
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It  is  of  interest  to  analyze  and  compare  the  term  containing 
the  field  integration  and  j  ((o)  separately.  Curves  for 
I(<W).  Hoo),   and   ^P(£  )  ^(Co,e  )(L€      are  given  in  Figure  9 
for  the  full  Coulomb  interaction  and  in  Figure  10  for  the 
dipole  approximation.  As  is  demonstrated  in  these  figures, 
the  term  j(^),  which  is  equal  to  four  times  the  K,,  matrix 
element,  is  importemt  only  in  the  line  center.  Its  region 
of  importance  is  slightly  greater  in  the  dipole  approximation 
than  in  the  full  Coulomb  interaction  case.  This  difference 
is  primarily  responsible  for  the  shoulder  obtained  from  the 
full  Coulomb  computations  and  illustrated  in  Figure  8. 

A  treatment  of  the  line  shape  problem  that  uses  a  full 
multi-pole  expansion  for  the  electron-atom  interaction  to 
all  orders  in  the  coupling  constant  has  been  carried  out  by 
M.  E.  Bacon,  K.  Y.  Shen,  and  J.  Cooper, ^^'^^  in  the  case  of 
the  Lyman-alpha  transition,  all  multi-poles  higher  than 
quadrupolar  terms  vanish  due  to  the  orthonormality  of  the 
atomic  states. 

The  major  features  of  the  calcxaation  by  Bacon  et  al. 
can  be  summarized  as  follows t 

1.  The  S  matrix  is  summed  to  all  orders  in  the  coupling 
constant . 

2.  Time  ordering  is  included  in  the  S  matrix. 

3.  The  perturbing  electron  trajectories  are  represented 
by  classical  paths. 

4.  Minimum  and  maximum  impact  parameter  cut-offs  are 
used  to  account  for  strong  collisions  and  electron  correlations 


Figure  9.  Functions  KCn)),  j(u)),  ^  (tt*,g)  for  the  full 
Coulomb  interaction  case* 
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respectively • 

5,     A  minimum  velocity  cut-off  is  employed. 

A  comparison  of  the  present  theoretical  restilts  with 
those  of  Bacon  et  al •  are  presented  in  Figure  11  for  a 

temperature  of  20,000*  K  and  an  electron  density  of  1.0  x 

17    3 
10  /cm  .  Only  the  center  and  near  ving  of  the  profile  are 

plotted.  The  agreement  is  fairly  good  for  wavelengths 

o 

greater  than  0.4  A.  However,  the  profiles  differ  consider- 
ably for  A)<   <,  0.4  A. 

The  full  Coulomb  treatment  presented  in  this  disser- 
tation uses  a  strong  collision  cut-off.  However,  reasonable 
variations  in  the  value  of  the  strong  collision  cut-off  does 
not  improve  the  agreement  between  the  two  theories  near  the 
line  center.  Also,  if  electron  correlations  are  accovinted 
for  by  a  cut-off  rather  than  by  the  use  of  the  electron  pair 
distribution  fvinction,  the  agreement  is  not  improved. 

In  the  quantum  mechanical  treatment  prerented  here,  the 
atom-perturbing  electron  potential  is  transformed  into  a 
momentum  representation  in  which  convergence  at  zero  momentum 
is  assured  by  other  momentum  dependent  factors  in  the  inte- 
grand for  ACt)  not  equal  to  zero.  For  zero  frequency  separa- 
tion the  manner  in  which  electron  correlations  are  treated 
assures  convergence  in  the  line  center  (see  Sections  II. 6  and 
II. 7  for  further  discussion).  The  major  convergence  factor, 
xP(q),  appears  explicitly  in  the  derivation  after  the  time 
integration  appearing  in  ^  (CO)  is  performed.   In  second 
order  theories  of  spectral  line  broadening  in  plasmas  time 


Figure  11.  A  comparison  of  the  Lyman-alpha  profile  using 
the  full  Coulomb  Interact: 
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ordering  can  be  significant.    This  may  mean  that  the  time 
ordered  integrations  in  higher  order  terms  do  not  provide 
convergence  factors  of  the  same  form  as  ^(q).  At  the 
present  time»  the  effects  of  time  ordering  in  higher  order 
terms  are  not  fxilly  understood. 

One  other  point  regarding  higher  order  terms  should  be 
mentioned.   In  the  present  treatment  the  W22  matrix  element 
appearing  in  ^(60)  is  approximately  1  percent  as  large  as 
the  W33  term.   In  the  fiai  CoiUomb  treatment,  the  states  \^ 
and  13")  ,  defined  in  Chapter  III,  are  not  mixed.  Hovever, 
they  will  be  mixed  if  higher  order  terms  in  the  perturbation 
expansion  are  retained.  Thus,  although  a  higher  order  term 
that  would  constitute  a  relatively  small  effect  on  the  second 
order  w^^  term,  may  through  mixing  of  states  in  higher  orders, 
have  a  relatively  large  effect  on  the  much  smaller  second 
order  Wjj  term. 

Another  effect  of  using  the  full  Coulomb  interaction  is 
to  reduce  the  asymmetry  in  the  line  profile  near  the  line 
center  as  is  illustrated  in  Figure  12  (compare  curves  1  and 
3).  This  behavior  is  a  consequence  of  the  differences  in 
magnitude  of  the  unshif ted  components  in  the  two  treatments. 
The  unshifted  components  in  this  second  order  theory  are 
broadened  only  by  the  perturbing  electrons.  It  can  be  shown 
that  the  electrons  are  mainly  responsible  for  the  asymmetry 
encoiintered  in  this  version  of  the  relaxation  theory t   a 
profile  computation  in  which  ^(CO)   is  symmetric  results 
in  a  symmetric  line  shape.  In  the  region  of  the  profile 
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where  the  asymmetry  is  large,  Jk6^/JjiT  ^,002   or  ^>  <.4  A, 

the  major  contribution  to  the  intensity  in  the  dipole 
approximation  arises  from  the  function  j((tf  )   given  in  equa- 
tion (5.3).  However,  in  the  full  Coulomb  treatment  the 
percent  contribution  from  the  function  j(6J)  to  the  intensity 
is  not  as  large  in  this  region.  Hence,  asymmetries  in  the 
function  J^£)  tend  to  be  amplified  in  the  dipole  approxi- 
mation compared  to  the  f\ai  Coulomb  treatment,  thus  resulting 
in  a  larger  asymmetry. 

A  comparison  of  a  line  profile  calculated  using  the 
full  Coulomb  interaction  with  the  experimental  data  and  with 
other  theories  appears  in  Figures  13  and  14.  The  agreement 
bet^-reen  the  profiles  appears  to  be  quite  good  over  the 
region  for  which  the  experimental  data  is  accurate  i   for  a 
wavelength  separation  greater  than  approximately  0.8  A.  The 
accuracy  of  the  experimental  data  is  limited  for  wavelength 
separations  less  than  0.8  A  due  to  self  absorption  effects. 
In  reference  6  it  is  estimated  that  in  the  experimental  data 
presented  here,  the  uncertainty  resulting  from  self  absorption 
corrections  can  be  as  large  as  10-15  percent  at  0.5  A  and 
even  greater  at  smaller  wavelengths.  Hence,  a  more  accurate 
procedure  for  performing  self  absorption  corrections  would  be 
necessary  in  order  to  determine  whether  the  shoulder  encoun- 
tered in  the  computations  presented  here  is  physically 
meaningful  or  T<hether  it  is  a  weakness  in  the  theoretical 
treatment.  If  the  latter  proves  to  be  the  case,  a  more 
extensive  analysis  including  higher  order  terms  in  the 
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perturbation  series  In  both  the  Ion  and  electron  Interaction 

terms  may  be  necessary. 
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(V.5)  ION-FIELD  EFFECTS  IN  THE  ELECTRON  COLLISION  OPERATOR 

The  Stark  splitting  of  the  energy  levels  is  directly 
proportional  to  the  ion-field  strength!  hence,  by  the 
analysis  presented  in  Section  V.2,  its  effect  should  be  most 
important  in  the  line  wings.  This  proves  to  be  the  case  as 
can  be  seen  in  Figure  12,  which  contains  plots  of  the  red 
to  blue  intensity  ratios.  Tmen  the  ion-field  is  included 
in  the  electron  collision  operator,  ^(<0),   the  asymmetry 
is  decreased  (compare  curves  1  and  4  in  Figure  12). 
Analytically  the  asymmetry  can  be  understood  from  an  analy- 
sis of  the  operator,  ^iCD).     In  the  wings  of  the  line  the 
major  cause  of  the  asymmetry  in  the  profile  arises  from  the 
factor  exp(XAfii>>^r-    ^mu'^^^^/^V^   appearing  in  Appendix  C. 
For  any  given  frequency,  larger  values  of  £  will  reduce  the 
relative  effects  of  asymmetry  in  this  factor.  Hence,  the 
asymmetry  in  the  line  profile  is  also  reduced. 

Physically  the  asymmetry  in  ^^{(0)   can  be  related  to 
the  nature  of  the  perturbing  electron-atom  interaction.  The 
following  heuristic  explanation  is  proposed.  The  energy 
difference  between  a  lower  energy  Stark  shifted  level  of  an 
excited  state  and  the  ground  state  energy  is  less  than  the 
energy  difference  between  the  unshifted  Stark  level  and  the 
ground  state.  Hence,  mixing  of  the  lower  Stark  state  with 
the  ground  state  is  increased,  thus  increasing  the  probability 


103 
of  a  transition*  On  the  other  hand,  the  higher  energy  Stark 
level  has  a  larger  spatial  extension  than  the  unshifted  level » 
thus  increasing  the  strength  of  the  perturbing  electron* 
atom  interaction  and  therefore  also  increasing  the  transi- 
tion probability.  These  two  effects  tend  to  counterbalance 
one  another »  but  are  not  expected  to  be  of  the  same 
magnitude  for  all  frequency  separations. 
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(V.6)  EFFECTS  OF  LOWER  STATE  INTERACTIONS 

In  the  dlpole  approximation  the  effect  of  electron 
perturbations on  the  lower  or  final  state  is  identically 
zero  since  the  matrix  element  of  the  atomic  dipole  operator 
is  zero  for  the  final  state.  Hovever,  in  the  full  Coulomb 
interaction  treatment*  the  lower  state  matrix  element  of 
the  perturbing  electron-atom  interaction  is  not  zero.  The 
effect  of  including  lower  state  interactions  is  illustrated 
in  Figure  15.  In  this  particular  curve  no  strong  collision 
cut-off  is  employed.  When  a  cut-off  is  employed  the  magni- 
tude of  this  effect  is  only  slightly  reduced.  Inclusion  of 
the  lower  state  interactions  leads  to  a  narrowing  of  the  line 
profile.  This  indicates  that  perturbations  that  decrease 
(increase)  the  energy  of  the  initial  level  also  decrease 
(increase)  the  energy  of  the  final  level*  thus  narrowing 
the  spread  in  energy  differences  between  initial  and  final 
states  resulting  in  a  narrowing  of  the  line  profile. 
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(V.7)  THE  EFFECTIVE  DISTRIBUTION  FUNCTION  RESULTS 

The  effect  of  electron  symmetry  on  the  Lyman-alpha 
profile  vas  Included  through  the  use  of  an  effective  distri- 
bution function,  p(x)t  for  the  perturbing  electrons  (see 
Section  II •2)*  The  computations  vere  restricted  to  the 
dlpole  approximation  for  computational  convenience* 

The  contribution  to  the  profile  from  matrix  elements 
W33  and  W44  appearing  In  ^  {60)     Is  dominant  for  AJ.  >1.0  A. 
Since  these  tvo  elements  are  of  the  same  fomit  ve  discuss 
the  effects  of  electron  symmetry  only  for  the  w,,  matrix 
element.  A  comparison  of  this  term,  calculated  using  an 
effective  distribution  function*  vlth  the  same  term  calcu- 
lated using  various  strong  collision  cut-offs  Is  given  In 
Figure  7.   By  plotting  W33,  at  several  frequencies,  as  a 
function  of  strong  collision  cut-off  and  comparing  these 
results  vlth  those  obtained  using  p(x),  an  effective  strong 
collision  cut-off  was  obtained.  The  resulting  effective 

cut-off  lies  bettireen  14  and  20  Bohr  radii,  which  Is  larger 

1  2  22 
than  cut-offs  currently  being  used.  *  *    Computations  using 

this  effective  cut-off  agree  with  those  using  p(x). 

The  present  calculations  using  p(x)  Increase  the 

Intensity  In  the  line  center  and  decrease-  the  Intensity  In 

the  line  wings.  The  amount  of  chzmge  Is  approximately  a 

9  percent  decrease  at  0.6  A  and  an  8  percent  Increase  at  3  A. 
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A  comparison  vith  other  treatments  is  best  accomplished  by 
plotting  the  ratio*  Ri 

where  I  ((•))  Is  the  Holtsmark  profile.  Such  comparisons 
appear  In  Figure  16.  For  wavelengths  greater  than  1,5  A, 
the  ratio  obtained  using  p(x)  does  not  Increase  as  rapidly 
as  In  other  treatments,  thus  resulting  In  better  qualitative 
agreement  with  the  experlmentea  data  for  1.5  A  than  achieved 
by  other  treatments.  It  Is  expected  that  the  agreement  may 
Improve  If  higher  order  terms  In  the  perturbation  expansion 
for  ^{^)   are  retained. 

The  red  to  blue  wing  ratio  exhibits  a  small  decrease 
In  magnitude  as  Is  Illustrated  In  Figure  18.  A  decrease  In 
the  ratio  Is  expected  since  the  asymmetric  electron  collision 
term  Is  smaller  In  magnitude  when  the  effective  distribution 
function  Is  used)  hence,  the  relative  effects  of  asymmetries 
In  this  term  are  correspondingly  reduced. 


4J 

« 
^ 

«H 

fl) 

IM 

B 

« 

•H 

h 

0) 

A 

a 

■P 

s 

01 

0) 

c 

0) 

•H 

.c 

S 

■p 

0 

•o 

« 

s 

4J 

I 

& 

% 

8 

u 

« 

■M 

01 

H 

•0 

•H 

d) 

IH 

•H 

0 

*M 

v« 

•H 

a 

§ 

1 

0) 

a 

J3 

;3 

4J 

1 

0 

g 

+> 

I 

g 

a 

-H 

4J 

id 

U 

IM 

§ 

0 

IM 

c 

fi 

o 

0 

n 

•H 

•H 
»4 

t 

« 

A 

a 

•H 

E 

M           • 

0 

•P         10 

u 

n      4J 

•H      Id 

< 

•0     -d 

• 

\0 

H 

0) 

^ 

01 

h 


110 


I 

•        1 

n 
■P 

1 

0 

5 

1 

0) 

+» 

Id 

CO 
0) 

M 

■g 

■P 
Id 

H 

^•^ 

•H 

•0 

•H 

X 

•h 

•h 

>-» 

•H 

H 

0 

p« 

§ 

Id 

u 

4> 

a 

1 

C 

Id 

g 

1 

0) 

1 

11  o 

U 

•H 

c 

M 

l*o 

10 

H 
•0 

•H 
H 

^ 

— 

V  "5^ 

U  Ob 

1 

1 

Id 

0) 

\oo 

h) 

Q 

OT 

o 

^  \         o 

• 

t 

O 

1 

•o 

o\\ 

°  V 

01 
•H 

>,     /vv^ 

5^    / 

oNSL 

s   / 

"^ 

*  / 

^ 

V 

4i       ' 

\ 

c—* 

^ 

0  X 

. 

•H"-* 

^^ 

\v 

0)  n 

s 

K   0 

o 
o 


u 

I 


J-^ 


(C»7)"l/((*?)I 


Ill 


(V.8)  STATIC  WING  RESULTS 

In  the  derivation  of  the  Holtsmark  mlcrofield,  ^^.i^)! 
lon-lon  interactions  and  the  effects  of  electron  screening 
are  neglected*  For  large  values  of  the  microfield,  i^ich 
are  important  in  the  line  vings*  Pi.(^)  provides  an  adequate 
description  of  the  microfield  function*  To  support  this 
claim*  reference  is  made  to  more  exact  treatments  of  the 
microfield  by  Baranger  and  Mozer,   and  Hooper.   ' 
For  large  values  of  the  field*  their  results  reduce  to  the 
Holtsmark  distribution*  Thus*  one  %rould  expect  the  ratio 
of  the  more  exact  theoretical  profiles  to  the  profile 
computed  in  the  static  limit  using  the  Holtsmark  distribu- 
tion should  approach  a  constant  for  large  frequency  separa- 
tions* For  the  static  wing  formula  developed  in  Chapter  IV 
the  asymptotic  value  of  the  ratio*  R  defined  in  Section 
(V*7)«  is  chosen  to  be  two*  If  the  dynamic  electron  terms 
do  not  have  the  proper  static  limit*  the  asymptotic  limit 
of  R  for  large  frequency  will  not  in  general  approach  this 
value*  Such  behavior  is  encountered  in  relaxation  theory 
computations* 

The  theoretical  results  for  the  static  wings  are 
compared  in  Figure  17*  If  one  assumes  that  the  experimental 
data  is  valid*  the  full  unified  theory  and  the  relaxation 
theory  in  the  dipole  approximation  appear  to  underestimate 
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the  value  of  R   for  small  vavelengths,  .^^  <  3  A.  This 
result  may  indicate  that  the  broadening  due  to  static  elec- 
tron fields  Is  underestimated  In  these  form\aatlons. 

In  performing  the  microf leld  convolutions,  better  agree- 
ment vlth  experimental  data  is  obtained  vhen  low-frequency 
microf lelds  are  used  for  both  the  electron  and  ion  micro- 
fields  than  when  either  high-frequency  or  HoltsmarK 
microf lelds  are  employed.  The  low-frequency  microf lelds 
account  for  electron  and  ion  correlations.  It  is  expected 
that  the  ion-ion  interactions  are  fully  shielded  (in  the 
Debye  sense)  by  the  elebtrons  in  the  plasma  whereas  the 
electron-electron  interaction  will  not  be  diielded  to  the 
same  degree.  Hence »  the  lo*r-frequency  convolution  should 
form  a  lower  bound  and  the  Holtsraark  convolution  an  upper 
bound  for  the  value  of  R   in  the  static  wing  limit. 
In  the  vicinity  of  the  electron  plasma  frequency, 

e 

^■^  <  2  A,  inclusion  of  higher  order  terms  in  the  time 
dependent  perturbation  expansion  given  in  equation  (4.12) 
is  expected  to  Improve  the  agreement  between  the  static 
wing  theory  and  the  experimental  data.  These  time  dependent 
terms  introduce  electron  dynamics  into  the  wing  formula. 
As  em  alternative  to  using  the  time  expansion,  one 


could  develop  a  theory  that  employs  time  dependent  mlcro- 

45 
fields   to  take  into  account  electron  dynamics.  However 

such  a  procedure  has  been  shown  to  be  computationally 

45 
difficult   since  it  is  necessary  to  evaluate  multi- 
dimensional Integrals  numerically.  Hence,  using  the 
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perturbation  In  time  developed  in  Chapter  IV  may  result  in 
a  more  tractable  form  for  computing  line  profiles* 

The  asymptotic  static  wing  formula  obtained  in  the 
unified  theory  contains  a  Holtsmarlc  microfield  distribution 
function  that  attempts  to  include  electron  correlations 
through  a  spatial  cut-off.  To  explore  the  validity  of  such 
a  treatment t  a  convolution  of  Holtsmarlc  microfields  computed 
using  a  spatial  cut-<off  similar  to  that  used  in  the  unified 
theory  was  performed*  The  results  are  almost  identical  to 
the  convolution  of  two  Holtsmarlc  microfields  without  a 
cut-off.  Hence*  it  is  concluded  that  the  convolution  of 
a  Holtsmark  distribution  with  a  cut-off  Holtsmark  distribu- 
tion will  be  insensitive  to  the  value  of  the  cut-off. 
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(V.  9)  THE  LYMAN-ALPHA  PROFILE  FOR  IONIZED  HELIUM 

In  order  to  illustrate  the  application  of  the  present 
theory  to  another  hydrogenic  system*  computations  vere 
performed  for  the  Lyman-alpha  transition  in  ionized  helium 

at  a  temperature  of  46(500*  K  and  an  electron  density  of 

17   3 
3*6  X  10  /cm  •  The  full  Coulomb  and  dlpole  treatment 

results  are  compared  in  Figure  18.  The  qualitative  features 

for  the  full  Coulomb  treatment  are  very  similar  to  the  Lyman- 

alpha  profiles  in  hydrogen,  that  is*  the  shoulder  is  fairly 

pronounced  and  the  agreement  betveen  the  two  theories  good 

in  the  line  vings*  Hovever,  it  is  of  interest  to  note  that 

for  this  particular  case,  the  shoulder  is  fairly  pronovinced 

in  the  dipole  treatment  also. 

At  the  present  time,  highly  accurate  experimental  data 

is  not  available  for  the  Lyman-alpha  transition  in  ionized 

46 
helium.    The  primary  problem  in  obtaining  accurate 

experimental  data  is  the  difficulty  encountered  in  measuring 

ultraviolet  radiation.  In  the  case  of  the  above  transition, 

the  vavelengths  involved  are  less  than  350  A. 
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(V.IO)  COMMENTS 

At  the  pre^ont  time  there  are  several  treatments  of  the 
line  broadening  problem  In  plasmas.  The  differences  In  these 
treatments  arise  mainly  from  the  different  methods  used  In 
evaluating  the  electron  collision  term*  Some  of  the  methods 
used  In  describing  the  electron-atom  Interaction  In  the 
electron  collision  term  arei 

a)  The  Interaction  between  the  perturbing  electron  and 
radiating  atom  Is  treated  as  a  dlpole  Interaction  and  a 
quantum  mechanical  perturbation  expansion  Is  developed  to 
second  order  In  the  coupling  constant.  »^2,23 

b)  The  Interaction  Is  described  by  a  dlpole  potential 
and  the  perturbation  expansion  Is  carried  out  to  all  orders 
In  the  coupling  constant.  Classical  paths  are  used  to 
describe  the  trajectories  of  the  perturbing  electrons. ^»^^ 

c)  A  complete  multi-pole  expansion  Is  used  for  the 
Interaction  auid  a  perturbation  expansion  Is  carried  out  to 
all  orders  In  the  coupling  constant  using  classical  paths  to 
describe  the  electron  trajectories. ^^'^^ 

The  vork  presented  In  this  dissertation  has  Introduced 
the  following  new  features! 

1)  The  Interaction  la  treated  by  an  exact  Coulomb 
potential  and  a  quantum  mechanical  second  order  perturbation 
expansion  employed. 
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2)  A  quasi-static  wing  formula  is  developed  that  is 

based  on  a  convolution  of  static  electron  amd  ion  microfields* 

3)  The  approach  listed  in  (a)  above  is  extended  to 
include  the  effects  of  electron  symmetry  in  an  approximate 
manner. 

The  differences  between  the  full  Coulomb  interaction 
treatment  and  the  other  treatments  listed  in  items  (a)  through 
(c)  above*  indicate  that  in  the  line  center  higher  order 
terms  in  the  coupling  constant  perturbation  expansion  can 
become  significant*  Hence*  it  would  be  interesting  to 
develop  a  theoretical  treatment  that  uses  a  full  Coulomb 
interaction  and  includes  all  orders  in  the  coupling  constant 
perturbation  expansion* 

Since  electron  symmetry  effects  are  shown  to  have  some 
effect  on  the  profile*  a  more  exact  treatment  of  the  effec- 
tive distribution  function  for  the  electrons  may  be  warranted* 
Such  a  treatment  should  take  into  account  the  fact  that  the 
radiating  atom  is  initially  in  an  excited  state.  It  is 
expected  that  the  probability  of  finding  an  electron  a 
distance  x  from  the  atom  will  depend  on  the  particidar 
excited  state  of  the  atom.  It  may  still  be  possible  to 
take  into  account  the  effects  of  electron  symmetry  in  an 
average  way  through  an  effective  cut-off. 

Another  problem  that  shottld  be  considered  in  developing 
a  static  wing  formula  is  the  computation  of  a  mlcrofield  for 
a  two-component  system  in  which  the  two  species  have  charges 
of  different  sign.  It  is  expected  that  the  convolution  used 
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in  this  presentation  Is  adequate  In  the  vings  of  the  profile, 
/^>  S  5*0  A.  HoT^ver,  for  smaller  vavelength  separations, 
a  microfield  distribution  calculated  with  a  more  exact 
treatment  of  the  electron-ion  interaction  may  be  required. 

Yet  another  extension  of  the  present  work  vould  be  to 
reduce  the  higher  order  terms  in  the  time  dependent  pertur- 
bation expansion,  presented  in  Chapter  IV,  into  a  tractable 
form  emd  to  compute  line  profiles  including  these  higher 
order  terms*  It  is  hoped  that  carrying  out  this  procedure 
will  lead  to  a  ving  profile  that  vill  overlap  treatments  of 
the  line  center. 

In  conclusion,  there  are  still  many  unresolved  questions 
regarding  the  line  broadening  problem.   I  hope  that  this 
vork  vill  add  to  the  understanding  of  the  line  broadening 
phenomenon  in  plasmas. 


APPENDICES 


APPENDIX  A 
THE  EFFECTIVE  POTENTIAL  V 


s 


In  this  appendix  an  effective  potential  Including 
exchange  Is  developed  for  a  perturbing  electron-atom  system* 
The  electron-proton  potential  Is  obtained  from  the  radial 
distribution  function  g^pC  f  »J^)»  defined  by, 

(A.l) 

where  j>  s  1/AT,  V^p  Is  the  effective  potential  or  potential 
of  mean  force,  /"  «  (2fP/4;'m).  vhere  ^  Is  the  reduced 
mass.  The  effective  Interaction  operator  Oj((<«))  Introduced 
In  Chapter  II  represents  a  free  electron  scattering  term. 
Hence,  the  summation  In  equation  (A.l)  vlll  be  restricted 
to  the  scattered  states  only.  These  states  are  described 
by  Coulomb  vave  functions.   '    Thus,  equation  (A.l)  becomes 

o 
(b   =  Bohr  radius  (A.  2) 


F       Is  the  Coulomb  vave  function. 
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The  above  expression  has  been  computed  by  Barker. ^^  The 
integration  over  jf^   is  trtancated  at  a  sufficiently  high  value 
to  assure  four  place  accuracy  in  determining  V*.,  The 
suRunatlon  over  ^  is  trtincated  -when  F^   becomes  less  than 
1.0  X  10""   .  For  ^a  0  the  Coulomb  functions  are  generated 
using  methods  derived  by  Frobert,^^  The  higher  order  Cou- 
lomb wave  functions  are  generated  by  formulae  given  in 

reference  47.  As  seen  from  equation  A.l»  the  potential  V 

ep 

is  the  logarithm  of  the  radial  distribution  function. 

The  effective  interaction  between  electrons  is  obtained 
in  a  similar  manner.  However*  for  the  spin  «  1  case,  the 
sum  over  J^  is  restricted  to  odd  values  of  J^  .  All  values 
of  jf  are  allowed  for  the  spin  «  0  case.  The  effective 
electron-electron  potential  is  obtained  by 

V^e  -   f^   <<l/4)g„p.^^^+  <3/4)9«p-up)         (A.3) 

where  the  subscripts  indicate  the  spin  orientation.     The 
potential  V^  defined  in  Chapter  II  is  the  stim  of  V^^  plus 
the  electron-proton  potential  obtained  from  equation   (A.2). 

^-^ee*       -^"^^ep),  (A.4) 

Thus 9      P      becomes 

Pe  •  ^'^(-Vae) 

C  exp(-pAg)exp(-Vg)^  (A.5) 

The  effective  density  matrix  given  in  equation  (A.5)  is 
used  in  evaluating  ^4  (<0)      .  For  large  values  of  x 
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(x  greater  than  50  Bohr  radii),  exp(-Vg)  reduces  to  one  as 
expected.  A  non  linear  least  squares  fit  is  performed  on 
the  function  exp(-V^)  using  exponential  functions  allowing 
one  to  perform  the  spatial  integration  occurring  in  ^  (CO) 
analytically.  Since  the  integration  range  is  small,  the 
numerical  error  introduced  by  using  this  procedure  should 
be  negligibly  small.  It  is  estimated  that  the  error  due 
to  the  numerical  procedures  used  in  evaluating    ^ (a> ) 
is  less  than  1  percent. 


APPENDIX  B 
PROPERTIES  OF  THE  OPERATORS  Lq  AND  P 

Five  equations  involving  the  operators  Lq,   P  and  L.  are 
derived  in  this  appendix*  These  relations  are  required  in 
reducing  the  perturbation  expansion  for  %^(Ca>)  to  a  useful 
form* 

1«  Let  A  represent  a  general  operatoi^  then^ 

<'PA>  «  <A>^  (B.l) 

Proof  of  (B.l) I 

<PA>  «  Tr^  <PAf^'> 

«  Tr^  (  feTV^)^ 

Since  the  term  Tr^CA  p  )  does  not  depend  on  the  elec- 
tron coordinates  one  has 

<PA>  ■  Tr3(A  Pg)  «  0  Q.E.D. 

■vhere  the  normalization  condition  Tr^  ^e  "  ^  ^^^   been 
explicitly  used* 

2*  The  commutator  relationship*   fL^fPl,  is  examined 
next*  In  order  to  evaluate  this  commutator  the  matrix 
elements  of  Lq  and  P  are  required*  The  system  under 
consideration  is  a  composite  system  consisting  of  electrons 

126 


127 
and  atoms.  LBtj\>S    represent  atomic  states  and  •<^p 
represent  electron  states.  The  matrix  elements  of  p   are 

''here  ^  Is  the  Kronecker  delta  and 

Note  that  the  matrix  elements  of  P  are  written  in  terms  of 
a  diagonal  and  non  diagonal  part. 

Next  consider  the  matrix  elements  of  PA  betveen 
composite  states.  One  has 


%ysf  =  <'-^-"iJV^)/^P>=<^«<i,ji:<wiAiv>ivrp>^ 


(PA) 

(B,4) 
Using  (B.3)  one  can  show  that 

''^Up  - ^f^p4*^<'-iie,)2A,,^^.^    ,3.3, 

Consider  next  the  matrix  elements  of  L  where  L  is  the 
Liouville  operator.  The  matrix  elements  of  L  are  given  by 

jimm^^'m'  aham'   i^m'  Ma,'    At^jui' .  (b.6) 

Hence*  for  a  composite  atom~electron  system 

^-^^     mWm'^'       ^"*"P/  "^-^^i*  ^"''f'  '  (B.7) 
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Using  equations  (B.6)  and  (B,7)  the  matrix  elements  of 
the  compound  operator,  (PL^)  can  be  determined.  These  are 

(B.8) 
Lq  Is  defined  by  equation  (2.30)  to  be 

vhere  Hq  «  H^^  +  v^^ji^. 

Using  equations  (2,33),  (2.34),  (2.35)  and  (B.6),  one  has 

where  :ft"^^^  ~  ^  ""  w  , 

The  results  of  equations  (B.9)  and  (B.IO)  are  substituted 
Into  equation  (B.B).  The  sums  over  ^'^  *  m*  &*   can  be  easily 
performed  using  properties  of  the  delta  functions.  The 
results  are 

^  ^"      ^'(^'i      .    (B.ll) 

A  similar  procedure  Is  followed  to  determine  the  matrix 
elements  of  (LqPA) . 
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•^  '  (B.12) 

Equations  (B«ll)  and  (B.12)  are  used  to  obtain  the 
matrix  elements  for  the  commutator  T^qvP'!  * 

•*P    /Wtf/H-  OM^,  ^w  ^wV^  (B.13) 

If  one  assumes  F.a  Is  diagonal*  e.g.»  If  one  neglects 
the  back  reaction  In  the  electron  density  matrix »  then 

[Lq,P]  «  0,  (B.14) 

This  Is  not  a  good  approximation  for  strong  collisions. 
Hence*  care  must  be  taken  whenever  (B.14)  Is  used. 

3.  Using  the  commutlvlty  of  P  and  L  one  can  verify 
equations  (2.39)  and  (2.40) •   One  has 

(l-P)Lp^  «  Lod-P)  p^  +   (l-P)L^f?^ 

-   (1-P)  L,  p^  (B.15) 

Hence* 

<A(1-P)L>  «   <^(1-P)L^'^  ^ 
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Using  (B.l)  one  can  write 


<L(1-P)  A>  ■  <P  L(1«P)  A> 

«   <PL^(1-P)  a;>   +  <LqP(1-P)  A>  . 


But  P(l-P)  ■  P  -  P^  ■  P  -  P  n  0|  hence. 


<L(1-P)  A)    -  <L^(1-P)  A>  .  (B.16) 


4»     The  next  relationship  to  be  proved  is 
l^iW-  CI-P)L^]''  s  [><W-L^]"(|-pJ  (B.17) 


Expanding  the  left  hand  side  in  powers  of  L.,  one  has 
or,  since   [Lq,  pI  «  0, 


where  the  idempotent  property  of  P  has  been  used.  Hence, 
Pi^-  (l-P)^,]  =  O^-L,"]  (I-P)^  (B,18) 


5.   In  a  similar  manner  as  above  the  following  identity 
can  be  provedi 

<B.19) 


APPENDIX  C 
A  COMPUTATIONAL  FORM  FOR  C04  (^  )^  .   i  i 

(C.l)  THE  FULL  COULOMB  TREATMENT 

A  computational  form  for  the  operator  ^JiiuJ)"^ 
occurring  in  equation  (2.54)  is  derived  in  this  appendix. 
The  derivation  is  carried  out  for  both  the  dipole  approxi- 
mation and  the  ftdl  CoulonA)  interaction  between  the  radiating 
atom  and  the  perturbing  electron* 

In  equation  (2.54)  the  operator  ^^^^  )/J^y  V>f  '    ^" 
given  in  terms  of  four  summations.  The  summands  of  the  four 
terms  are  related  to  one  another  through  permutation  of  the 
indices.  Hence,  it  will  be  sufficient  to  consider  only  one 
of  these  terms  in  detail.  In  the  derivation  that  follovs, 
the  third  term  in  equation  (2.54)  will  be  considered  and 
be  referred  to  as  F((^).  Using  equation  (2.49)  for  Jl(A<ti) 
one  has 


(C.l) 


where  in  this  appendix  STZ     ,  is  defined  to  be 
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The  eigenvalues  of  the  operator  H^  when  operating  on  the 
wave  functions   <^ )  and  <p|   are  Jf^U)      and  J\  U)^ 
respectively,  with  these  eigenvalues  the  exponential, 
exp(-  t  (m)      X    )»  can  be  vrltten  In  operator  notation.  Hence 
equation  (C.l)  becomes 


FMc  -/  i/^(^^^,,.o  I  «\^^^  Ip)^ 


>p 


.* 


IC  .  2  / 


T-ihere  V^^  ,  represents  matrix  elements  of  V^^^  between  the 
atomic  states  x\  and  ^ . 

The  Inverse  Fourier  transform  of  function  F(4;)  has  been 
treated  In  reference  23  using  the  dlpole  approximation  for 
V^^  and  an  uncorrelated  single  electron  distribution  func- 
tion defined  as  lA»  where  V  Is  the  volume  of  the  system. 
This  derivation  and  that  presented  In  reference  23  are 
essentially  the  same  up  to  equation  (A.19)  In  that  reference. 
Using  equation  (A.19)  In  reference  23  and  equation  (C.2),  one 
obtains 

(C.3) 
where  a  »  4i^tti  /-^T,  s  «  t^^T^  ,  and  q  represents 
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moment  von  transfer  for  the  perturbing  electron* 

The  s  integration  can  be  performed  using  the  following 

t    «.    ,  47,50 
integral 6 I 

J  Si^{l^)j^i-c^)   =   WF^  D(l>A*rS-)         (C.5) 


vhere  D(x)  is  Davson's  integral  defined  by 


The  s  integration  yields  a  function  of  qi 


»• 


Thus  performing  the  time  integration  the  expression  for  F(6i)), 


equation  (C*3),  becomes 


.-k  /-:»   4-^ 


(C.7) 
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The  y  integration  is  performed  next.  The  potential 

Vjjg(y)  can  be  vrltten  as 

^  •^      X  *  (C.8) 

vhere  j  Is  the  perturbing  electron  Index.  Furthermore,  ve 
may  write 


i<»/>  =  -'^ [iip^\  ~  \\,.\) 


vhere  R  and  yj  are  the  position  operators  for  the  atomic 
and  perturbing  electron  respectively,  and  the  nucleus  Is 
assumed  to  be  located  at  the  origin.  The  Integration  over 
"y .  Is  just  the  Fourier  transform  of  the  Coulomb  potential 
vhlch  Is 

For  the  total  potential  defined  In  equation  (C.8)  the  "y" 
Integration  yields 

^^  (C.IO) 
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Inserting  these  results  Into  the  expression  for  F{^),   one 
finds 

tJ 
X  <v|  Z  (x^(-/*..rt/*)-d)Aa  IV';). 


(C.ll) 

There  are  tvo  distinct  terms  in  the  above  expression 
for  F(a>),  j  »  k  and  j  p'  K.  The  j  ■  k  term^  labeled  Fq,  is 
considered  first.  The  subscripts  are  suppressed  for  compact- 
ness. 

Since  all  the  perturbers  under  consideration  are 
identical*  the  summation  over  j  can  be  replaced  by  N,  the 
total  number  of  electrons*  times  the  j  ■  1  terms 


<v  \^(-^'l'  RM)-0/5*  \  vr'>  T  M^;)A'A^^^. 


(C.12) 
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All  the  integrations  over  qj  except  i  ■  l  can  be 
readily  performed  using  the  properties  of  the  Dlrac  delta 
function.  Hence f 

><   ^^ ''*•"'»-•  (C.13) 

we  define 


(C.14) 


where  p  ( x,)  Is  a  reduced  distribution  function  and  Is 
Interpreted  as  representing  the  probability  of  finding  a 
perturbing  electron  at  x   regardless  of  the  location  of 
the  other  N-1  electrons.  An  approximation  for  p(x:)  that 
takes  Into  account  the  effects  of  electron  symmetry  Is 
developed  In  Appendix  A. 

The  potential  (l/j^  -t)     -  '/i4,\    )    occurring  In 
equation  (C.13)  can  be  vrltten  In  terms  of  an  Inverse 
Fourier  transform 


(C.15) 


Using  the  definition  of  p(x^)   in  equation  (C.14)  and 


equation  (C.15),  Pq  becomes 


137 

(C.16) 
vhere  subscripts  on  x  and  q  have  been  suppressed  for 
compactness  In  vrltlng*  and  n^   Is  the  electron  density » 

The  atomic  matrix  elements  of  the  operator 
(expdq'R/^  )>1)  are  evaluated  In  Appendix.  D*  Let 

fvw-'  -  <^ '  -^^  '*^'  ^/^ )  -  M \r' >.  (C.17) 

In  terms  of  G  .  i  ,  F^  becomes 
Vir    0 

The  angular  part  of  the  x  Integration  Is  performed  next, 
choosing  tk.^   along  the  q»-q  direction* 

)(.  S,v(^/|l-^'l/if)/^,^.|i|/^^  X^^gfA^.  (C.19) 
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The  angular  integrations  involving  q*  and  q  are 
performed  next.  The  coordinate  system  for  the  q*  integra- 
tion is  chosen  so  that  q^*   lies  along  the  q  direction.  After 
carrying  out  the  angular  integration  in  q',  the  angular 
integration  in  q  is  trivial.  F.  becomes 

o 

(C.20) 

For  sufficiently  large  x,  x  S  x  <5^  50  Bohr  radii, 

p(x)<b:l*  Thus  it  is  convenient  from  a  computational  stand 

point  to  split  the  x  integration  up  into  tvo  regions,  x  /  x 

m 

and  X  V  X  ,  where  p  -Ji  1  for  x  S  x  .  Hence 


-^ 


-h  1  (f»'<'*-)~  ^^)  S'^*A/(^^M)Sc*A/(y.Vy^)^^^^ 


21) 


The  first  integral  is  carried  out  by  vriting  the  sine 
terms  in  exponential  form  and  using  the  integral  definitions 
of  the  Dirac  delta  functional 
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A, 


(C.22) 

By  substituting  equation  (C«22)  into  equation  (C.20) 
and  performing  the  q*  integration,  one  can  show  that 


+  PTERM  (C.23) 

where  PTERM  is  defined  as 


PTEBM     „  ^ 


6 

(C.24) 

After  performing  a  minor  chzuige  in  variables*  the  first 
term  in  expression  (C.23)  is  evaluated  nvunerically. 

The  second  term,  which  contains  p(x) ,  has  been 
examined  by  using  a  non  linear  least  squares  fit  for  p(x)  • 
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The  X  integration  can  then  be  performed  analytically.  The 

remaining  tivo-dimensional  q»  q'  integral  contains  oscillatory 

functions  and  is  poorly  behaved,  hence  it  is  very  difficult 

to  evaluate.  However »  in  the  dipole  approximation  the  q* 

integration  can  be  performed  analytically*  Thus,  the  detailed 

effect  of  using  the  effective  distribution  function  is 

examined  only  in  the  dipole  approximation*  For  the  full 

Coulomb  interaction  treatment »  only  an  order  of  magnitude 

estimate  is  obtained  for  the  second  term*  and  this  by  means 

of  an  effective  cut-off  (see  Section  V.5). 

The  term  j  y^  k*  labeled  F^,  occurring  in  equation  (C.ll) 

is  evaluated  next.  Electron  correlation  effects  will  enter 

into  the  theory  via  this  term. 

Due  to  the  identical  nature  of  the  perturbing  electrons 

N(N>1)  identical  terms  will  result  T^en  the  summations  are 

carried  out.  Making  the  approximation  N(N-l)  »  N  and 


setting  j  ■  1,  k  <*  2  one  has 


(C.25) 


Defining  p(x^,X2)  by 


■f<^.,^^^  =  v*i'--$PK-'^-W^»--'''^*' 


(C.26) 
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and  using  the  Inverse  transform  of  V   one  can  show  that 


vhere  G  ^p'   and  G^/^'  are  defined  in  equation  (C.17). 

The  angular  Integrations  In  q"and  q*  are  performed  next, 
q^  and  q^»  are  taken  to  lie  along  the  x  and  "Sfj  directions 
respectively*  One  obtains 

It   correlations  between  electrons  are  neglected*  p(x  ,x  KbS'O 
and  F^  Is  zero  due  to^  the  angular  Integrations  over  x,  and 

Xj. 

In  terms  of  the  pair  distribution  fianctlon,        p(x  .x-) 
may  be  written  as 


(C.29) 
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and 


The  spatial  integrations  are  carried  out  next.  A 
transformation  to  relative  coordinates  is  made  and  the 
definitions  appearing  in  equation  (C.30)  are  used  to  obtain, 

(C.31) 
The  angular  integration  in  y  yields  a  factor, 

iS-/^)^  ^^  '         <C.32) 

Substituting  equation  (C. 32)  into  equation  (C.31)  and 
performing  the  angular  integration  in  x  yields, 

<;/!^w>^  ^/»  -^^-^  ^Ao  s-^»^/^) 

^  (C.33) 

vhich  is  further  reduced  by  performing  the  x  integration 
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(C.34) 


S 


i^ere  O  (x)  Is  the  Dlrac  delta  function*  When  these  results 
are  Inserted  Into  F^  and  the  q*  Integration  carried  out  one 
obtains  ^  -f-        J^ 


(C.35) 

Combining  equations  (C«35)  and  (C.23)  one  obtains  P  for 
the  case  of  full  Coulomb  interaction  between  the  radiating 
atoms  2md  perturbing  electrons » 


tMjfjTt*  r   ^     ^  f^ti^A.' 


-s  ■       --  %^ 


1l'     -     i+«»1/ii)*  J    , 


PTERM  ,  (C,36) 


The  above  expression  is  evaluated  by  numerical 
techniques. 
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C.2  THE  DIPOLE  APPROXIMATION 


If  the  dlpole  approximation  is  used,  the  potentials 
V(x)  and  V(y)  in  equation  (C.7)  are  replaced  by  the  dlpole 
interaction  between  the  atom  and  perturbing  electron.  Then. 

^^  -Apr-  I^^3^>  S  ?,(L.-^'A;^  W^'^i'M 


a/ 


X 


Mi 

(C.37) 

Equation  (C.37)  also  yields  ttfo  terms,  1  =  j,  and  1  f'  j. 
The  1  «  j  term,  labeled  Pj^^  in  the  following  derivation,  is 
evaluated  first.  Following  a  procedure  similar  to  that  used 
in  evaluating  F  one  obtains 

0  . 

(C.38) 
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The  y  Integration  is  performed  using  the  coordinate 
system  Illustrated  belov. 


The  results  of  this  integration  are 


A     n^^» 


F   - 


[\i(f-)[j^-/'^^)'i.M'^/^' 


X  ^^X#•'?/^  )  C.«(  p^^^^  )/^  ^Jj^J^ 


(C.39) 


The  X  integration  can  be  performed  by  a  similar  proce- 
dure yielding  the  following  expression  for  F  i 


(C.40) 
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The  angular  part  of  the  q  Integration  is  performed  uaing 
the  coordinate  system  Illustrated  belov. 


'  W  %\/* 


This  integration  yields 

(C.41) 

A  similar  procedure  is  folloved  for  the  1  /  j  term. 
The  details  are  rather  lengthy  but  straightforward  and  will 
be  omitted.  The  restating  expression  for  the  correlation 
term  (see  Section  C.l),  Fj^,,  is 
iJ^M^  /?  ,  .1?  , 


be 


5^flT 


(C.42) 
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Combining  equations  (C.42)  and  (C.4l)  one  has,  in  the 
dipole  approximation, 

^  =  ~^^'  ^4  5r'^^  <*/^> 


(C.43) 


If  one  follows  the  commonly  accepted  procedure  of  setting 
p(x)  equal  to  unity  and  using  a  strong  collision  cut-off,  X* 


the  expression  for  F  _  becomes 


(C.44) 


^^%^w/^ 


F    =  -   "*"*        v(r 


3>a"r 


^Jr^)   [^'W2P/./i)/(^>^./;t:) 


'•/('.^a^/>3^^]    ^^, 


(C.45) 

Computer  programs  have  been  written  for  the  various  ap- 
proximate forms  of  F  given  in  equations   (C.36),    (C.43)   and   '^C.45). 


APPENDIX  D 
MATRIX  ELEMENTS  OF  EXPd^^^r)  -  1 

The  atomic  matrix  elements  of  the  operator  exp(i A  t)-!, 
required  In  the  evaluation  of  //k{io)^  ,    ,    ,   are 

evaluated  In  this  appendix. 

The  parabolic  vave  functions  for  a  hydrogen  atom  are 
given  by 

$  <'<.>.♦) = r  '^  <-—'-'-'>  ^ — __f 

^A*^  L  C/»/»4^')'/v»7r((sy./>M,0/(44-s-/);) 

(/vmI 


(D.l) 

vhere  aQ  Is  the  Bohr  radius,  L^  are  the  associated  Laguerre 

47 
polynomials,   and.  In  this  Appendix,  n,  s,  m  are  the  para- 
bolic quantum  numbers* 

The  transformation  equations  between  the  parabolic  and 
spherical  coordinate  systems  are. 
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X-  4T^  <**C<P)  (D.2a) 

y  -  >|Acr  S,'AfC<i>)  (D.2b) 

z  -  0-6t)/a.  (D.2c) 

r  «  O+U)/  a.  (D.2d) 

Ji^    -(C>-f-^)/f.)  i>^H^4>  .                 (D.2e) 

If  a  coordinate  system  is  chosen  so  that  the  z  compo- 
nent  of  "r   lies  along  the  ^  direction » 

exp(iA'r)^^  exp(i^z)-^  exp(i  A(  A -H  )/2)^   (D.3) 

In  the  following  derivation  ^^   and  V  vill  be  used  to 
represent  the  atomic  states.  Using  this  notation  and  the 
above  expression  for  exp(i^*r)  one  has 


«    o 


xa^^«v-h)  ^(-if  <ir^^^)  L  ^'^''"-'-^ 


l-vKxc 


|>i*«^  \  \M^  \ 


(D.4) 


ISO 
where 


(D.5) 
The  integration  over  p     yields  a  factor  2  IT  &jy^    /mi 

^  "^   IT 

vhere  cp^  .  Is  the  Kronecker  delta*  The  remaining  integra- 
tions depend  only  on  the  absolute  value  of  m.  Therefore, 
for  convenience  (ml  win  be  denoted  as  m  where  m  is  assumed 
positive. 

A  transform  of  variables  is  performed  in  equation  (D«4)t 

"  (D.6a) 

Then 

X  C<~^  vy  liJi;^^)  r^J]^^)  f^A.^)  ^^^^ 


s 

(D.7) 
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vhere 


/M^f  ^^^ 


(D.8) 


In  the  no-quenching  approximation  one  requires  matrix 
elements  for  vhich  /M^  »  /^^m   n,  where  n  is  the  principal 
quantum  number  of  the  initial  state,  in  this  case  Jl    - 
-^\r"  ^»  ■"**  ^  "   ^»o"/2.  Hence, 


<D.9) 
vhere 
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Equation  (D.9)  can  be  written  in  terms  of  the  folloving 
integral 

(D.IO) 
resulting  in  the  following  expression* 

'  -*  •       (D.ll) 


The  integral  h(i»m,j»k)  can  be  reduced  to  a  closed  form 
that  can  be  easily  programmed  for  machine  computation*  First 

note  that  the  generating  function  for  the  associated  Laguerre 

38 
polynomials  is  ^tf 

rCM  +  r-t-O-^    H   /  JVM  - 


T  P(Tf/)  ^ 


{•(^"■--i 
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Carrying  out  the  indicated  derivatives  In  the  above  expres- 
sion yields  the  following  svunmatloni 


(D.12) 


tfhere    r^  j    Is  the  generalized  binomial  coefficient  defined 
by 

\q  I  ql(p-q)l     ' 

Inserting  the  above  results  Into  equation   (D.IO)   leads  to 


(T-zy^-  ^*"'^(-^> 


(D.13) 


The  X  Integration  can  be  performed  using  the  Integral 
derived  from  formulae  appearing  In  reference  50, 

(D.14) 
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One  finally  arrives  at  the  following  closed  expression  for 

h(ifm»  j»1c).  j^  I 


J>=o 


(n)-^  2c-.f(r.)^, 


^i  p.. 


j-t-r 


c  A^ 


(9.15) 


A  computer  program  has  been  written  to  evaluate  the 
above  expression i  because  it  is  diagonal  in  the  magnetic 
quantum  number*  m,  the  matrix  exp(i-A*r)-lt  vill  at  most 
be  of  order  n. 


APPENDIX  E 
USEFUL  MATRIX  RELATIONSHIPS 

Tvo  matrix  relationships  useful  in  line  shape  computa- 
tions are  derived  in  this  appendix. 

A  complex  matrix,  M,  can  be  written  as,  M  ■  A  +  iB, 
where  A  and  B  are  real  matrices.  Its  inverse  is  written  as, 
M~  B  C  +  iD.  Expressions  for  C  and  D  can  be  derived  in 
terms  of  A  and  B  and  their  inverses.  By  definition, 

(A  +  IB)  (C  +  ID)  «  1 

where  1  is  the  identity  matrix.  The  above  expression  yields 
tiro  matrix  equations  when  multiplied  out 

AC  -  BD  B  1 
BC  -  AD  «  0. 

Solve  for  C  emd  D. 

C  «  (A  +  BA"*^B)"^ 
D  at  -(B  +  AB"-^A)"^  . 

Assume  that  A  can  be  written  as 

A  «  A*  +  otE 

where  E  is  diagonal  and  <9^  is  a  scalar  parameter.  Then, 
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(B  +  AB"^A)    «    (B  +    (A«   +      AE)b'^(A»   +      AE) ) 


B  +   A*b"^A«   +  «<(AEB"^A'   +  A'b'^AE)    - 


+  ol^aeb'^ae 

«  Bl   +    o(B2  +     ot?B3 

vhere  the  definitions  of  Bl  and  B3  follow  from  the  derivation. 

When  the  above  relationships  are  applied  to  the  line 
shape  problem*  many  of  the  matrix  products  can  be  performed 
outside  of  the  ion-field  integration  routine,   thus  saving  a 
considerable  amoiint  of  computer      time* 


APPENDIX  F 
DIPOLE  MATRIX  ELEMENTS  IN  THE  PARABOLIC  REPRESENTATION 

The  dlpole  matrix  elements  of  the  atomic  dlpole  operator 
are  evaluated  directly  In  the  parabolic  representation  In 
this  appendix.  The  derivation  Is  partially  based  on  material 
presented  In  Chapter  II  of  reference  18. 

The  z  component  of  the  dlpole  operator  Is  considered 
first.  In  the  parabolic  coordinate  system  one  has 

2  «  (  >  -  M,  )/2 

vhere  y  and  C\  are  the  parabolic  coordinates  as  defined  In 
Appendix  O.  Following  a  similar  procedure  as  that  used  to 
compute  the  matrix  elements  of  exp(c  b(  V  -^))  in  Appendix 
D,  one  can  derive  the  following  expression  for  the  matrix 
elements  of  z  i 


fr    o 


^      V  IT 


s  s;^    -^  A_         Mi    ^ 


where  the  terms  Involved  In  this  expression  are  defined  In 

Appendix  D. 
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Defining  the  following  Integral 


m  *'*-^ 


one  can  vrlte 


<M\i\\ 


"3M.*»r 


The  matrix  elements  of  x  and  y  components  can  be 
obtained  In  terms  of  an  operator,  fjT  t  defined  as 


After  performing  the  p   Integration  (see  equation  D.9), 
one  obtains I 


M,     V 


/H  <f  -^Z 


vhere 
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and 


fi.  p^!  ("-S^- 1-8,1-0 (  S^{  (M-^-  l«i,l  -■)(  ' 


The  Integral  I  is  evaluated  in  reference  18. 
Using  the  above  results,  one  arrives  at  the  follovlng 
matrix  elements  for  the  Lyman-alpha  transitions  In  hydrogen i 


(x) 


<^. 


0  C    0    0    c 

C  0  -1,5  -1,5    0 

0  -1.5   0     0  -1.5 

0  -1.5    0    0  -1.5 

C  0-1.5  -1.5   0 


(F.l) 


(y) 


OJ^ 


■  0  C    0    0  c 

-C  0   1.5   1,5  0 

0  -1.5    0     0  1.5 

0  -1.5    0    0  1.5 

-C  0  -1.5  -1.5  0 


(F.2) 


(2) 


1    0 

0 

c 

-c 

0 

0 

0 

0 

0 

0 

c 

0 

-3 

0 

0 

-c 

0 

0 

3 

0 

0 

0 

0 

0 

0 

(F.3) 

vhere   o^  is  the  Bohr  radius,  c  «  .5267482  and  the  state 
labels  are 
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1  B      |100> 

2  B      |20-1> 
3=       |210> 
4S      |200> 
5S      |200>      , 

For  the  Lyman-alpha  transitions  the  above  matrix  elements 
could  have  been  obtained  directly  from  the  wave  functions 
without  developing  a  general  expression  for  the  matrix 
elements.  However,  such  a  procedure  would  be  undxily  compli- 
cated for  the  higher  transitions. 

To  obtain  the  matrix  elements  for  helium  ion  transi- 
tions* one  replaces  the  Bohr  radius  by  Ct^/2, 


BIBLIOGRAPHY 

1.  H.  R.  Griem,  Plasma  Spectroacopv  (McGraw-Hill  Book  Co., 
New  York,  1964). 

2.  E.  w.  Smith,  C.  R.  Vidal,  and  J.  Cooper,  J.  Res.  Nat. 
Bur.  Stand.   73A,  389  (1969). 

3.  M.  Baranger,  Atomic  and  Molecular  Procesaes.  D.  R,  Bates, 
Ed.  (Academic  Press,  Inc.,  New  York,  1962). 

4.  R.  Zwanzig,  Lectures  in  Theoretical  Physics,  w.  E.  Downs 
and  J.  Down,  Eds.  (Interscience  Publishers,  Inc.,  New 
York,  1961) I  J.  Chem.  Phys.  33,  1338  (i960) i  Phys.  Rev, 
124,  938  (1961) I  Physica  30,  1109  (1964). 

5.  R.  Kubo,  Lectures  in  Theoretical  Physics,  W.  E.  Downs 
Etfid  J.  Down,  Eds.  (Interscience  Publishers,  Inc.,  New 
York,  1959). 

6.  R.  C.  Elton  and  H.  R.  Griem,  Phys.  Rev.  135,  A1550  (l964)i 

7.  G.  Boldt  and  W.  S.  Cooper,  Z.  Naturforsch  19A,  968  (1964). 

8 .  H .  Bremermann ,  Distributions,  Complex  Variables,  and 
Fourier  Transforms  tAddlson>Wesiev  !{»ubiiahing  ffn..  Tnn.^ 
Reading,  Massachusetts,  1965). 

9.  E.  w.  Smith,  Dissertation,  University  of  Florida  (1966). 

10.  L.  Galatry,  Phys.  Rev.  122,  1218  (1961). 

11.  E.  w.  Smith,  J.  Cooper,  w.  R.  Chappell,  and  T.  Dillon, 
to  be  published. 

12.  E.  W.  Smith  and  C.  F.  Hooper,  Jr.,  Phys.  Rev.  157,  126 
(1967).  

13.  E.  w.  Smith,  J.  Cooper,  and  C.  R.  Vidal,  Phys.  Rev.  185, 
140  (1969).  

14.  C.  F.  Hooper,  Jr.,  Phys.  Rev.  165,  215  (1968) i  Phys. 
Rev.  169,  193  (1968). 

15.  C.  F.  Hooper,  Jr.,  Phys.  Rev.  149,  77  (1966). 

16.  A.  C.  Kolb  and  H.  R.  Griem,  Phys.  Rev.  111.  514  (1956). 

161 


162 

17.  H.  R.  Grlera,  A.  C.  Kolb,  and  K.  Y.  Shen,  Phys.  Rev.  116. 
4  (1959).  

18.  P.  Kepple,  Report  #831,  University  of  Maryland  (1969). 

19.  H.  A.  Lorentz,  Proc.  Acad.  Sci.  Amsterdam  8,  591  (i906). 

20.  V.  weisskopf,  Z.  Physics  75,  287  (1932). 

21.  M.  Levis,  Phys.  Rev.  121,  501  (l961). 

22.  H.  R.  Griem,  Phys.  Rev.  140,  A1140  (1965). 

23.  E.  W.  Smith,  Phys.  Rev.  166,  102  (1968). 

24.  A.  A*  Barker,  Aust.  J.  Phys.  21,  121  (1968). 

25.  F.  w.  Sears,  Thermodynamics,  the  Kinetic  Theory  of  Gases, 
and  Statistical  Mechanics  (Addison-wealev  Puhi '1  «i>i]ng  rrs., 
Inc.,  Reading,  Massachusetts,  1964). 

26.  u.  Fano,  Phys.  Rev.  131,  259  (1963)  i  Rev.  Mod.  Phys.  29, 
74  (1957) I  Lectures  on  the  Many-Body  Problem.  E.  R. 
Caianiello,  Ed.  (Academic  Press,  Inc.,  New  York,  1964) . 

27.  N.  Davidson,  Statistical  Mechanics  (McGraw-Hill  Book 
Co.,  Inc.,  New  York,  1962).       ' 

28.  K.  Huang,  Statistical  Mechanics  (John  Wiley  and  Sons, 
Inc.,  New  York,  1965). 

29.  A.  A.  Barker,  J.  Chem.  Phys.  55,  1751  (1971). 

30.  A.  Ohno,  Linear  Stark  Broadening  by  Charged  Particles 
(Quantum  Chemistry  Group,  Uppsala  U.,  Uppsala,  Sweden, 
1963),  Preprint  No.  109. 

31.  B.  Mozer,  Dissertation,  Carnegie  Institute  of  Technology 
(1960).  a-y 

32.  E.  Merzbacher,  Quantum  Mechanics  (John  Wiley  and  Sons, 
Inc.,  New  York,  1960). 

33.  S.  S.  Schweber,  An  Introduction  to  Relativistic  Quanttmi 
Ffeld  Theory  (Harper  and  Row,  Publishers,  Inc.,  New  York, 
1961) . 

34.  E.  w.  Smith,  Private  Communication,  Fall  1971.. 

35.  E.  W.  Smith,  Phys.  Rev.  Letters  18,  990  (1967). 

36.  J.  T.  O'Brien,  Dissertation,  University  of  riorida  (1970). 

37.  L.  A.  Jones,  J.  R.  Grieg,  T.  Oda,  and  H.  R.  Griem,  Phys. 
Rev.  4,  833  (l97l). 


163 

38.  P.  M.  Morse  and  H.  Feshback,  Methods  of  Theoretical 
Physics  (McGrair-Hlll  Book  Co.,  Inc.,  New  York,  1953). 

39.  J.  w.  Dufty,  Phys.  Rev.  187,  305  (1969). 

40.  R.  M.  Wilcox,  J.  Math.  Phys.  8,  962  (1967). 

41.  M.  Bacon,  F.  K.  Shen  and  J.  Cooper,  Phys.  Rev.  188,  50 
(1969). 

42.  M.  E.  Bacon  and  D.  F.  Edvards,  Colorado  State  Univer- 
sity, Physics  Department  Technical  Report  No.  TR-68-1 
(unpublished),  1967. 

43.  J*  T.  Godfrey,  C.  R.  Vidal,  E.  w.  Smith  and  J.  Cooper, 
Monograph  No.  121,  (slat.  Bur.  of  Stand.,  Washington, 

D.  C.  1971). 

44.  M.  Baranger  and  B.  Mozer,  Phys.  Rev.  115.  521  (1959). 

45.  L*  J.  Roszman,  Dissertation,  University  of  Florida  (1971). 

46.  J.  R.  Grieg,  H.  R.  Griem,  L.  A.  Jones,  and  T.  Oda,  Phys. 
Rev.  Letters  24,  3  (1970) i  Phys.  Rev.  A4,  833  (l97l). 

47.  M.  Abramovitz  and  I.  A.  Stegun,  Handbook  of  Mathematical 
Functions  (U.  S.  Government  Printing  Office,  Washington, 
D.  C,  1964). 

48.  N.  F.  Mott  and  H.  S.  w.  Massey,  The  Theory  of  Atomic 
Collisions  (Oxford  Univ.  Press,  London  and  New  York, 
1949). 

49.  C.  A.  Frobert,  Rev.  Mod.  Phys.  27,  399  (1955). 

50.  S.  M.  Selby  and  B.  Girling,  Standard  Mathematical  Tables 
(The  Chemical  Rubber  Co.,  Cleveland,  Ohio,  1965). 

51.  T.  L.  Hill,  Statistical  Mechanics  (McGrav-Hill  Book  Co., 
Inc.,  New  York,  1960). 


BIOGRAPHICAL  SKETCH 

Joseph  Edvard  Whalen  vas  born  July  7,  1944,  In  Boston, 
Massachusetts.  In  June,  1962,  he  was  graduated  first  In 
his  class  from  Belllngham  High  School,  Belllngham,  Massachu- 
setts.  In  June,  1966,  he  received  the  degree  of  Bachelor 
of  Science  with  a  major  In  physics  from  Worcester  Polytechnic 
Institute  In  Worcester,  Massachusetts.  He  entered  the 
Graduate  School  of  the  University  of  Florida  In  September, 
1966.  He  worked  as  a  graduate  teaching  assistant  In  the 
Department  of  Physics  and  Astronomy  until  January,  1971. 
In  January,  1971,  he  took  on  the  duties  of  Scientific 
Programmer  for  the  Department  of  Physics  and  Astronomy  at 
the  University  of  Florida,  a  post  which  he  has  held  up  to 
the  present  time.  Throughout  his  enrollment  at  the  University 
of  Florida  he  has  pursued  his  work  toward  the  degree  of 
Doctor  of  Philosophy. 

Joseph  Edward  vjhalen  Is  married  to  the  former  Judith 
Anne  Campbell  and  Is  the  father  of  tt^o  children.  He  Is  a 
talented  commercial  and  jazz  pianist  and  has  led  and 
participated  In  many  musical  groups  since  his  early  high 
school  days.  He  Is  a  member  of  Sigma  Xl  and  the  American 
Physical  Society. 


164 


I  certify  that  I  have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality,  as 
a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 


ialrman 
Associate  Prd'fessor  Af  Physics 


I  certify  that  I  have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  ftilly  adequate,  in  scope  and  quality, 
as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 


<2.c 


ilSr 


R.  C.  I8l( 

Associate  Professor  of  Physics 


I  certify  that  I  have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 


E.  s.-T^ 


-*vv.><,^v 


B.  S*  Thomas 

Assistant  Professor  of  Physics 


I  certify  that  I  have  read  this  study  and  that  in  my 
opinion  it  conforms  to  acceptable  standards  of  scholarly 
presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 


J.tff.  Dufty     J   5 
AsJplstcmt  Frofestfbr  of 


Physies 


I  certify  that  I  have  read  this  study  and  that  in  my 
opinion  it  conf ozms  to  acceptable  standards  of  scholarly 
presentation  and  is  ftilly  adequate,  in  scope  and  quality, 
as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 


C.  p.  Luehr 

Assistant  Professor  of  Mathematics 


This  dissertation  was  submitted  to  the  Department  of 
Physics  and  Astronomy  in  the  College  of  Arts  and  Sciences 
and  to  the  Graduate  Council*  and  was  accepted  as  partial 
fulfillment  of  the  requirements  for  the  degree  of  Doctor  of 
Philosophy, 

Dean,  GTraduate  School  — — — — — — 


u 

■liaffii 


